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ABSTRACT 

Although  parts  of  this  paper  may  be  applicable  to  various 
other  location  finding  systems,  radio  direction  finding  was 
selected  for  development  because  of  the  extensive  literature 
available  on  the  subject.   Using  a  modification  of  the  navi- 
gation coordinate  system,  equations  in  spherical  trigonometry 
are  developed  by  analogy  with  graphical  solution  procedures. 
Ellipticity  of  the  earth  is  corrected  for  by  use  of  geocen- 
tric latitude  in  all  calculations.   The  least  squares  solu- 
tion is  found  by  a  combination  of  second-order  gradient 
searches  and  Golden  Section  Searches  on  the  corrected  sphere. 
A  Chi-Square  contour  confidence  region  is  presented  with 
defense  against  the  empty  region  arguments  and  with  verifi- 
cation of  its  reliability  of  target  containment.   Sensitivity 
of  location  error  to  network  layout,  number  of  measurements 
available  and  intercept  geometry  factors  is  examined  by 
computer  simulation,  demonstrating  the  predictive  value  of 
these  statistics.   It  is  shown  that  the  correlation  between 
location  error  and  the  size  of  a  reliable  confidence  region 
is  too  minimal  to  have  predictive  value. 


TABLE  OF  CONTENTS 

I.  INTRODUCTION ■ 7 

II.  MATHEMATICAL  DESCRIPTION  OF  SPATIAL 
RELATIONSHIPS  9 

A.  GLOBAL  COORDINATE  SYSTEM  9 

B.  GEOCENTRIC  LATITUDE  10 

C.  THE  POLAR  TRIANGLE 14 

1.  Point  Azimuth  Solution  15 

2.  Point  Azimuth  and  Signal  Azimuth 
Difference 17 

III.  THE  FIX  SOLUTION 21 

A.  OBJECTIVE  FUNCTION  22 

B.  SECOND-ORDER  GRADIENT  SEARCH  22 

C.  GOLDEN  SECTION  SEARCH  27 

D.  FORTRAN  FIX  PROGRAM 31 

1.  Program  Description ■ 32 

2.  Results 34 

IV.  TARGET  LOCATION  CONFIDENCE  REGION  35 

A.  SIGNAL  AZIMUTH  FRROR  DISTRIBUTIONS  35 

B.  THE  OBJECTIVE  FUNCTION  AS  A  CHI-SQUARE 
STATISTIC 36 

C.  CONFIDENCE  LEVELS  AND  REGION  BOUNDARIES  37 

D  .   FORTRAN  CONFIDENCE  REGION  PROGRAM  47 

1.  Program  Description  48 

2.  Results 49 


V.  INTERCEPT  GEOMETRY  51 

A.  PLOTTING  POSITION  LINFS  IN  SOLUTION 

REGION 51 

B.  MEAN  ACUTE  ANGLE  OF  INTERSECTION 5  3 

C.  MEAN  ARC  STATISTIC 56 

VI.  SIMULATION  PROCEDURE  58 

A.   GEOCENTRIC  AND  GEOGRAPHICAL  LATITUDE  59 

VII.  LOCATION  ERROR  ANALYSIS  WITH  FIXED  DEGREES 

OF  FREEDOM 60 

A.  COMPARISON  OF  ERROR  ON  NATIONAL  AND 

GLOBAL  NETWORKS  60 

B.  ERROR  VERSUS  MEAN  ACUTE  ANGLE  OF 
INTERSECTION  62 

C.  ERROR  VERSUS  MEAN  ARC  STATISTIC 69 

VIII.  LOCATION  ERROR  ANALYSIS  WITH  VARYING  DEGREES 

OF  FREEDOM 76 

A.  NATIONAL  NETWORK  76 

B.  GLOBAL  NETWORK 80 

IX.  CONCLUSION 86 

COMPUTER  PROGRAM  88 

COMPUTER  OUTPUT  97 

LIST  OF  REFERENCES 108 

INITIAL  DISTRIBUTION  LIST  109 

FORM  DD  1473 111 


LIST  OF  FIGURES 

1.  Global  Coordinate  System  10 

2.  Relationship  Between  Geographical  and 

Geocentric  Latitude  11 

3.  Polar  Triangle  with  Signal  Azimuth  14 

4.  Azimuth  and  Vertex  Angle 18 

5.  Azimuth  Differences  21 

6.  Reducing  the  Search  Interval  28 

7.  Searching 34 

8.  Bounded  Signal  Azimuths  42 

9.  Bernoulli  Random  Frror  45 

10.    Right  Spherical  Triangle  52 

11a.   90  Degree  Acute  Angle  of  Intersection  55 

lib.   10  Degree  Acute  Angle  of  Intersection  55 

12.  Relating  Error  to  Emitter  Distance  56 

13.  Location  Error  Versus  Degrees  of  Freedom  61 

14.  National  Network  with  N  Equal  Three: 

Intersections  63 

15.  National  Network  with  N  Equal  Five: 

Intersections  64 

16.  National  Network  with  N  Fqual  Seven: 

Intersections ■ 65 

17.  Global  Network  with  N  Equal  Three: 

Intersections  66 

18.  Global  Network  with  N  Equal  Five: 

Intersections  67 

'19.    Global  Network  with  N  Equal  Seven: 

Intersections  68 

20.    National  Network  with  N  Equal  Three:   Arcs  70 


21.  National  Netv/ork  with  N  Equal  Five:   Arcs 71 

22.  National  Network  with  N  Equal  Seven:   Arcs  72 

23.  Global 'Network  with  N  Equal  Three:   Arcs 73 

24.  Global  Network  with  N  Equal  Five:   Arcs 74 

25.  Global  Network  with  N  Equal  Seven:   Arcs 75 

26.  National  Network:   Area  Versus  Error  77 

27.  National  Network:   Error  Versus  Degrees  of 

Freedom 78 

28.  National  Network:   Error  Versus  Intersections  79 

29.  National  Network:   Error  Versus  Arcs  81 

30.  Global  Network:   Area  Versus  Error  82 

31.  Global  Network:   Error  Versus  Degrees  of 

Freedom ■ 83 

32.  Global  Network:   Error   Versus  Intersections  84 

33.  Global  Network:   Error  Versus  Arcs 85 


I.   INTRODUCTION 

Various  procedures  have  been  described  for  estimating  an 
emitter  location  from  a  set  of  measured  signal  azimuths 
reported  to  a  control  station  by  dispersed  direction  finding 
(DF)  sites.   Such  procedures  range  upward  in  complexity  from 
visual  estimates  off  a  plot  of  great  circle  arcs  on  a 
gnomonic  projection  of  the  operating  area  of  potential 
targets.   Other  procedures  exist  for  describing  confidence 
criteria  about  the  location  estimate.   The  present  paper  was 
motivated  by  an  interest  in  confidence  regions  and  in  loca- 
tion error  analysis.   Such  studies  require  either  that  the 
procedure  used  be  free  of  approximating  assumptions,  or  that 
the  variability  of  the  approximations  be  included  in  the 
confidence  criterion.   The  former  course  was  preferred  for 
simplicity  and  consistency  of  estimates.   This  paper  presents 
an  emitter  location  estimate  procedure  based  only  upon  known 
DF  site  coordinates  and  signal  azimuths.   The  procedure 
described  in  sections  II  and  III  converges  quickly  to  a 
least  squares  solution  and  does  not  require  an  initialization 
point  estimate  of  the  location  from  external  sources  (Pope, 
1971) .   The  procedure  does  not  require  a  rejection  of  out- 
liers (McCalla,  1970)  in  order  to  converge  to  the  least 
squares  solution.   The  approach  is  free  from  assumptions  in 
the  categories  that  the  earth  is  flat  in  the  solution  region, 
that  the  earth  is  a  perfect  sphere,  that  great  circle  arcs 


are  well  enough  approximated  by  straight  lines  of  position, 
and  that  azimuth  errors  are  well  enough  approximated  by 
parallel  displacements  of  such  lines  of  position  (Daniels, 
1951) . 

Section  IV  describes  a  Chi-Square  contour  confidence 
region  procedure.   Factors  from  intercept  geometry  are 
discussed  in  section  V.   Section  VI  is  a  description  of 
procedures  used  to  simulate  the  problem  for  the  purpose  of 
the  location  error  analyses  described  in  sections  VII  and 
VIII.   An  alternative  to  confidence  regions  for  predicting 
location  accuracy  is  discussed  in  the  concluding  section. 


II.   MATHEMATICAL  DESCRIPTION  OF  SPATIAL  RELATIONSHIPS 

Locations  and  spatial  relationships  of  interest  in  direc- 
tion finding  and  target  fixing  can  be  expressed  completely 
in  spherical  trigonometry.   The  advantage  of  such  expression 
is  that  each  mathematical  step  has  its  simple  graphical 
counterpart.   This  paper  contains  notational  and  descriptive 
departures  from  other  direction  finding  models  based  on  spher- 
ical trigonometry.   Such  notational  changes  follow  those  used 
in  geodesy  for  describing  locations  and  spatial  relationships 
and  those  used  in  textbooks  on  spherical  trigonometry  for 
developing  equations  to  express  such  relationships.   Fquations 
are  restricted  to  sine  and  cosine  functions  only  in  order  to 
simplify  calculating  first  and  second  order  derivatives  of 
somewhat  lengthy  equations . 

A.   GLOBAL  COORDINATE  SYSTEM 

The  coordinate  system  is  a  slight  modification  of  the 
system  used  in  mapping.   Longitude,  denoted  A,  is  measured 
from  zero  at  the  prime  meridian  eastward  through  360  degrees. 
Latitude,  denoted  <j>,  is  measured  from  zero  at  the  equator 
northward  through  ninety  degrees  and  southward  through  minus 
ninety  degrees.   Azimuth,  denoted  8,  is  measured  from  zero 
at  true  north  clockwise  through  360  degrees.   Such  a  coor- 
dinate system  is  sufficient  to  represent  any  point  on  the 
globe  as  well  as  the  great  circles  which  represent  ideal  wave 
propagation  paths.   Mathematical  descriptions  of  spatial 


relationships  must  account  for  the  longitude  discontinuity 
at  the  prime  meridian  and  for  the  azimuth  discontinuity  at 
true  north.   A  great  circle  arc  with  azimuth  angle  indicated 
is  sketched  on  the  coordinate  system  in  Figure  1. 
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Global  Coordinate  System 
FIGURE  1 

B.   GEOCENTRIC  LATITUDE 

Geographical  latitude,  used  in  navigation  and  mapping, 
is  based  upon  astronomical  observations  and  mathematical 
deductions  that  do  not  correct  for  the  ellipticity  of  the 
earth.   Consequently,  geographical  latitude  is  slightly 
greater  in  magnitude  than  the  arc  on  a  meridian  from  the 
equator  to  a  location  with  vertex  angle  at  the  center  of  the 
earth.   Rules  of  spherical  trigonometry  can  be  more  accurately 
applied  to  describe  spatial  relationships  following  the  con- 
version from  geographical  latitude  to  geocentric,  or  earth- 
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centered,  latitude.  Assuming  that  waves  propagate  along 
great  circle  arcs,  then  no  information  is  lost  since  such 
paths  are  geometrically  independent  of  radii  to  the  earth 
center  from  both  the  signal  emitter  and  from  the  DF  site. 
The  maximum  discrepancy  between  the  two  types  of  latitude 
occurs  at  ±45  degrees  latitude  and  is  about  equal  to  11.5 
minutes  of  arc.  Although  the  conversion  equation  appears 
in  books  on  geodesy,  a  derivation  could  not  be  found,  so  the 
equation  was  derived  for  reference  purposes. 

In  Figure  2,  a  and  b  represent  the  semi-major  and  semi- 
minor  axes  of  the  ellipsoid.   The  point  P (x, ,y, )  represents 
any  point  on  the  ellipsoid.   The  line  OP  connects  P  to  the 
origin  at  0,  the  earth  center.   SP  is  perpendicular  to  the 


Equation  of  ellipse 


0 


Relationship  Between  Geographical  and  Geocentric  Latitude 

FIGURF  2 
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equatorial  plane,  the  x-axis .   NP  is  the  normal  line  to  the 
ellipsoid  at  P,  and  the  angle  PNS  is  equivalent  to  geographi- 
cal latitude.   Angle  PON  is  equivalent  to  geocentric  latitude, 
denoted  \p .      Ellipticity  of  the  quadrant  is  exaggerated  in 
order  to  distinguish  the  two  angles. 

From  Figure  2,  it  can  be  seen  that 

SP   =   OS    tan    \p ,    and    that   SP    =   NS    tan<|> , 

so    that 

tamj;  =  ||  tan<|>  (2-1) 

From  the  equation  of  the  ellipse, 

Y1  =   |  (a2  -  x2) '  (2-2) 

The  equation  of  the  normal  line  NP  is  given  by 

y  -  *i =  -  T^kr  (x "  xi}  (2_3) 


Differentiating  equation  (2-2)  to  solve  f  (x.) , 

f '  (Xj^)  =  -  -  x1  (a   -  xx) 

Substituting  in  equation  (2-3)  and  rearranging  terms, 

y  -  *i  =  bx^  (a2  -  xi}  2  (x  -  xi> 

2 

=   — X (x  -  x, )  ,  from  equation  (2-2) . 

b  x. 

Then  at  y  =  0,  y,  can  be  cancelled  from  both  sides,  and 
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x  -  x,  ,2 

1  _    b_ 

x,  2 

1  a 


From  Figure  2,  x,  =  OS ,  and  at  y  =  O,  x  =  ON,  so  that 


x,  -  x  =  NS. 

NS    b^ 

OS     2 
a 

Substituting  this  last  result  in  equation  (2-1) , 

h2 
tarnj*   =   ~  tancf>  (2-4) 

a 

The  conversion  is  usually  stated  as  a  function  of  earth 
eccentricity, 

1    /  2    K2^ 

e  =  (a   -  b  ) 


Then 


o     2    ,2 
2    a   -  b 
e   =  - 


2 

a 

K2 

=  1  -  — 

x     2 

a  . 

Therefore, 

b2        2 

— 2   -    1  -  e  ,  and  it  follows  that 

a 

tan^  =  (1  -  e2)  tan<J>  (2-5) 

Using  e  =  0.081813  from  the  International  Ellipsoid  of 
Reference  (reference  3) , 

tamjj  =  0.99327733  tancf)  (2-6) 
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C.   THE  POLAR  TRIANGLE 

The  mathematical  description  of  spatial  relationships  is 
based  upon  a  polar  triangle  constructed  from  the  great  circle 
arcs  (measured  by  the  angle  subtended  by  the  arc  at  the 
earth's  center)  connecting  a  DF  site,  any  other  point  and  the 
north  pole.   Such  a  triangle  is  sketched  in  Figure  3.   The  DF 
site  is  located  at  vertex  B.   Vertex  C  is  at  any  other  point. 
Vertex  A  is  at  the  north  pole.   The  arc  c  is  the  complement 
of  the  geocentric  latitude  of  the  DF  site.   Similarly,  the 
arc  b  is  the  complement  of  the  geocentric  latitude  of  the 
other  point.   Angle  A  is  equal  to  the  difference  in  longitude 
between  the  DF  site  and  the  point.   The  dotted  arc  included 
on  Figure  3  corresponds  to  a  signal  azimuth  with  angle  indi- 
cated by  6 . 

A 


Polar  Triangle  with  Signal  Azimuth 
FIGURE  3 
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Angle  A  corresponds  to  the  smaller  of  two  possible  polar 
triangles  connecting  any  two  points  and  the  north  pole. 
Visually,  there  is  no  problem  in  deciding  which  polar  tri- 
angle to  choose.    However,  when  this  triangle  contains  the 
prime  meridian,  the  discontinuity  at  longitude  360  =  longi- 
tude 0  results  in  assigning  the  value  of  the  vertex  of  the 
larger  polar  triangle  when  the  longitude  difference  is  taken 
The  mathematical  model  can  be  forced  to  correspond  to  the 
spatial  relationships  by  adding  the  proviso  that  if  the 
difference  in  degrees  of  longitude  is  greater  than  180 
degrees,  then  let  A  equal  360  degrees  minus  the  difference. 
This  restricts  the  size  of  angle  A  to  one  hemisphere.   The 
absolute  value  of  the  difference  can  be  taken  in  case  the  DF 
site  is  east  of  (greater  longitude  than)  the  point.   These 
relationships  are  summarized  below  with  subscript  s  desig- 
nating the  DF  site  and  subscript  t  designating  any  test 
point. 

c  =  90  -  \p  (2-7) 

b  =  90  -  ty  (2-8) 


A  = 


X.  -  A  I  ,  if  I A.  -  A    <    180 
t    s '      '  t    s '  — 


360  -  |xt  -  Aj,  if  |At  -  As|  >  180 


(2-9) 


1.   Point  Azimuth  Solution 

.  * 

The  value  of  arc  a  connecting  the  DF  site  and  test 
point  is  given  by  the  law  of  cosines  for  sides, 

cos  a  =  cos  b  cos  c  +  sin  b  sin  c  cos  A         (2-10) 
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There  is  no  quadrant  ambiguity  since  the  cosine  takes  nega- 
tive values  in  the  second  quadrant  and  arc  a  is  always  less 
than  or  equal  to  180  degrees  because  of  the  restriction  on 
angle  A  to  one  hemisphere.   Since  the  cosine  of  angle  A  is 
of  interest  in  determining  arc  a,  both  the  absolute  value 
sign  and  the  condition  can  be  removed  by  applying  a  reduction 
formula.   From  equation  (2-9) , 


cos  A  = 


cos   A .  -  A 
t     s 


cos  (360  -  I  A.  -  A   ),  if  I  A.  -  A  I  >  180 
1  t    s   '     '  t     s ' 


But,  through  reduction, 

cos  (360  -  I  A.  -  A  I)  =  cos  I  A ,  -  A 

1  t    s '         '  t     s ' 

The  absolute  value  sign  can  be  dropped  since  cos  (A   -  A  )  = 
cos  (A   -  A  ).   From  this  result  and  equation  (2-10) 

cos  a  = cos  b  cos  c  +  sin  b  sin  c  cos  (A.  -  A  )     (2-11) 

t    s 

The  vertex  angle  B  at  the  DF  site  can  be  found  by  applying 

one  of  the  laws  of  cosines  for  sides, 

cos  B  =  cos  b  -  cos  c  cos  a 
sin  c  sin  a 

Equation  (2-12)  is  valid  for  test  points  either  east  or  west 
of  the  DF  site  at  vertex  B. 

In  order  to  make  comparisons  with  a  signal  azimuth, 
angle  B  must  be  converted  to  an  azimuth  from  true  north, 
denoted  B.   Visually,  the  conversion  is  based  on  whether  the 
test  point  is  eastward  or  westward  of  the  DF  site  with  res- 
pect to  the  proper  p'olar  triangle.   Mathematically,  the 
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conversion  is  based  on  the  magnitude  of  the  longitudes  of 
the  DF  site  and  the  test  point,  but  must  also  account  for 
the  discontinuity  at  the  prime  meridian.   There  are  four 
general  cases  with  two  distinct  outcomes  as  illustrated  in 
Figure  4 . 

2 .   Point  Azimuth  and  Signal  Azimuth  Difference 

From  a  given  DF  site,  the  difference  between  the 
azimuth  of  any  test  point  and  a  signal  azimuth  is  the  smaller 
of  the  two  possible  outcomes,  depending  upon  whether  the 
difference  is  measured  clockwise  or  counter-clockwise.   The 
azimuth  difference  is  denoted  to,  and  the  azimuth  discontinuity 
at  true  north  requires  that  the  computation  be  conditional. 
With  a  clockwise  rotation  for  the  azimuths,  the  smaller 
angle  of  azimuth  difference  will  be  positive  if  3  is 
leading  and  negative  if  3  trails  0.   The  possible  outcomes 
and  their  conditions  are  as  follows. 

a.  If  | 3  -  8 |  <  180,  then  w  =  3  -  6 

b.  If  |  3  -  6|  >  180  and  3  >  6,  then  to  =  3  -  6  -  360 

c.  If  |3  -  6|  >  180  and  3  <  9,  then  to=3  -0  +  360 

(2-14) 
The  sets  of  conditions  on  computing  3 (B)  in  equation 
(2-13)  and  in  computing  to  (3,0)  in  equation  (2-14)  yield  six 
possible  outcomes  for  computing  azimuth  difference  as  a 
function  of  the  vertex  angle  B  and  signal  azimuth  0. 
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a.   A   <  A^_  and  A^_  -  A   <  180,  then  3  =  B 
st      t     s 


b.   A   >  X.  and  A  ^-A.  <  180,  then  3  =  360  -  B 
s     t      s     t 


c.   A   >  X^    and  A   -  A^  >  180,  then  3  =  B 
s   .  t      s     t 


d.   A   <  X.  and  X^    -    A   >  180,  then  3  =  360  -  B 
s     t       t     s 

Azimuth  and  Vertex  Angle 
FIGURE  4 
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u 


-< 


(  b  -  e 

B  -  360  -  e 

B  +  360  -  6 

-  B  +  360  -  9 

-  B  -  0 
V  -  B  +  720  -  0 

Two  utility  variables  may  be  introduced  in  order  to  express 
azimuth  difference  in  a  constant  form  with  values  dependent 
upon  the  existing  conditions. 


to  =  pB  +  q 


(2-15) 


(p,q)  =  < 


The  pair  of  utility  variables  can  take  on  any  one  of  the 
following  sets  of  values 

(    (1,0) 

(1,-360) 

(1,360) 

(-1,360) 

(-1,0) 
I  (-1,720) 

Concluding  development,  angle  B  can  be  expressed  as  a  differ- 
entiable  function  of  sines  and  cosines  of  known  values  of 
longitude  and  geocentric  latitude.   From  equations  (2-11) 
and  (2-12)  and  the  application  of  an  angle-difference  rela- 
tion, 
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B  =  cos 


(cos  b  -  cos  c  [cos  b  cos  c  +  sin  b  sin  c 

(cos  A^  cos  A   +  sin  A.  sin  A  )~|  ]  /  sin  c 
t      s        t      s  JJ   i 

sin  (cos   [cos  b  cos  c  +  sin  b  sin  c  (cos  A   cos  A 

v  V 

+  sin  A   sin  A  )]  J 

Substituting  the  values  for  arcs  b  and  c  from  equations  (2-7) 
and  (2-8)  and  using  the  reduction  formulas  cos  (90  -  x)  = 
sin  x  and  sin  (90  -  x)  =  cos  x, 


B  =  cos 


-1 


cos  ip 


(sin  ty,    -   sin  ty      [sin  iJj,  sin  ty      +    cos  ^, 

cos  ty       (cos  A.  cos  A   +  sin  A   sin  A  )] 

sin  (cos    [sin  ^   sin  ty      +  cos  ^   cos  iJj   (cos  A 
cos  A   +  sin  A   sin  A  )]  j 

By  multiplying  the  term  sin  ^   through  the  interior  brackets 

of  the  numerator,  writing  the  first  two  resulting  terms  as 

2  2 

sin  \/j   (1  -  sin  i>    )  ,    which  is  equal  to  sin  ip.  cos  ty    f 

cos  i>      can  be  cancelled  from  both  numerator  and  denominator, 
s 

B  =  cos     sin  i>,    cos  ty      -    sin  ip   cos  ty       (cos  A 

\  I  L-  o  o  l—  I— 


cos  A   +  sin  A,  sin  A  ) 
s        t      s 


sin   cos 


-1 


[sin  4>t 


sin  ty      +  cos  Uj,  cos  \p   (cos  A.  cos  A 
s        rt      rs        t      s 


+  sin  A   sin  A  ) 


(2-16) 
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Ill .   THE  FIX  SOLUTION 

In  Figure  5,  each  of  four  DF  sites  have  measured  azimuths 
of  a  signal  from  the  same  emitter  at  an  unknown  location. 
The  signal  azimuths  are  plotted  for  illustration  as  dotted 
arcs.   A  test  point,  P(X. ,ty.),    is  established,  and  azimuths 
from  each  DF  site  to  P  are  drawn  in  solid  arcs.   A  better 
fit  can  be  found  by  moving  P  somewhere  within  the  region- 
bounded  by  the  dotted  arcs . 

Given  a  number  n  of  DF  sites,  each  reporting  the  azimuth 
of  a  received  signal  with  respect  to  its  own  location  to  a 
net  control,  the  fix  problem  can  be  viewed  as  finding  the 


Azimuth  Differences 
FIGURE  5 
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point  on  the  globe  which  minimizes  the  sum  of  some  function 
of  the  azimuth  differences  co.,  for  i  =  1,2,..., n. 

A.  OBJECTIVE  FUNCTION 

Since  the  azimuth  differences  are  signed,  minimizing 
their  sum  may  produce  a  solution  that  lies  outside  the  region 
bounded  by  arcs  corresponding  to  signal  azimuths  when  large 
positive  and  large  negative  azimuth  differences  cancel  each 
other  out.   Minimizing  the  sum  of  the  absolute  values  of 
azimuth  differences  avoids  this  problem  of  positive  and 
negative  differences;  however,  since  the  partial  derivatives 
of  absolute  value  are  not  computationally  tractable,  solution 
techniques  are  limited.   Minimizing  the  sum  of  the  squared 
azimuth  differences  provides  an  objective  function  that  is 
dif f erentiable  and  reflects  absolute  magnitude.   Therefore, 
the  solution  of  the  fix  problem  is  taken  to  be  that  pair  of 
coordinates  (A  ,  ty    )    which  will 

n    2 
Minimize    E   w .  (A.  ,ip.)  (3-1) 

i=l   1    t   t 

B.  SECOND-ORDER  GRADIENT  SEARCH 

Dropping  the  t  subscript'  to  simplify  exposition,  the 
objective  function,  equation  (3-1),  can  be  represented  as 
f  (X  ,\p)  .      Given  an  initial  test  point  (A,ip),  a  truncated 
Taylor's  series  for  a  function  of  two  variables  can  be  used 
to  approximate  the  functional  values  of  other  test  points 
(A  +  AA,  \p   +  AtJ/)  . 
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f  (X    +    AX,    i>   +    Ai|;)    =    f  (X^)    +    A*    5f^^)    +    AX    8f^^} 


+  i(AX)2  s^iii  +  axa^  afgtjuti  +  i(A^)2  i!f%n 

2         8X2  8A3^      2         8^2 


Expressing  the  approximation  in  matrix  form  and  letting 
superscript  t  indicate  the  transpose, 

A,TfO\ 


'X    +    AX\  9X 

=    f  (A,!*;)11   +     (AA  ,Aip) 
4>   +   Aip  /  1    8f(A,i|i) 


8ip 


y 


/r_92f  (A, if;)       92f  (A,i|>)' 


+   |  (AX,  Ai|;) 


9X' 


8  A  3ip 


V 


2f(X,i|;)  2f  (A  ,^) 


AX^ 
Ai|;, 


9X9i|; 


3i|>' 


Letting  P.  denote  an  initial  test  point  (X,i|>)  ,  letting  R 

denote  the  vector  (AX,Ai[i)  ,  letting  V  denote  the  gradient 

2 

of  f  (P  . )  and  letting  V   denote  the  Hessian  matrix  of  f  at 

P  .  ,  the  approximation  can  be  expressed  as 

f  (P.  +  R)  =    f(P.)  +  RfcVf(P.)  +  -  RtV2f  (P.)  R 
3  J  3     2         3 

Then,  .. 

f  (P.  +  R)  -  f  (P.)  =  RfcVf  (P.)  +  ~   RfcV2f  (P.)  R        (3-2) 
3  3  3     2         3 

Starting  from  an  initial  test  point  at  P.,  the  search  pro- 
cedure can  be  understood  as  a  technique  to  find  a  vector  R 
such  that  f (P .  +  R)  is  the  minimum  of  the  objective  function. 
Such  a  vector  R  also  minimizes  the  left  hand  side  of  equation 
(3-2) .   Since  the  quadratic  form  on  the  right  hand  side  is  an 
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approximation  of  the  value  of  the  left  hand  side,  it  must 

be  near  a  minimum  value  when  the  left  hand  side  is  minimized. 

The  approximation  works  both  ways,  so  that  a  vector  R  which 

minimizes  the  right  hand  side  must  approximately  minimize 

the  left  hand  side.   The  right  hand  side  is  in  a  form  that 

yields  to  analysis  of  its  extreme  values  with  respect  to  R 

at  a  fixed  P . .   One  of  the  conditions  necessary  for  the 
3 

quadratic  form  to  be  at  a  minimum  is  that  its  gradient  be 
equal  to  zero.   Therefore,  it  is  necessary  to  find  a  vector 
R  which  fulfills  this  condition.   Taking  the  gradient  of 
equation  (3-2)  with  respect  to  R, 


R 


RfcVf  (P.)  +  i-  RtV2f  (P.)  R 

r       2  y 


=  v 


R 


RfcVf  (P.) 


+  V 


R 


|  RrV^f  (P  )  P. 


=  Vf  (P.)  +  \   V2f(P.)  R  +  \    RfcV2f  (P.) 
J     2       j       2         j 

=  Vf  (P  .)  +  V2f  (P  .)  R 
3  3 

Setting  the  gradient  of  the  quadratic  form  at  R  equal  to  zero, 

..  Vf  (P.)  +  V2f  (P  .)  R  =  0 
3  3 

Then, 


V2f  (P  .)  R  =  -Vf  (P.) 
3  3 

2 
And,  providing  that  V  f(P.)  is  invertible, 


R  =  - 


1   -1 


V"f  (P.) 


Vf  (P  ) 


(3-3) 
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R  specifies  both  a  direction  and  magnitude  in  which  to 
move  toward  an  extreme  value  of  the  objective  function. 
This  consequence  of  the  Taylor  approximation  can  be  improved 
by  applying  a  scalar  coefficient,  denoted  6,  to  R  in  order 
to  adjust  the  magnitude  and,  if  R  points  to  a  maximum  value, 
the  direction.   The  next  test  point  is  then  given  by 


P  ...  =  P  .  +  6R 
1  +  1     3 


Combining  equations  (3-3)  and  (3-4) , 


p\  i  =  P.  -  <5  [v2f  (P.jl   Vf  (P.) 
3+1  j  L     i'l  ■}' 


(3-4) 


(3-5) 


The  search  procedure  can  be  terminated  when  further 
improvements  are  of  little  consequence,  most  readily  deter- 
mined by  comparing  the  distance  between  locations  on  succes- 
sive iterations. 

In  order  to  solve  equation  (3-5) ,  it  is  necessary  to 
compute  the  inverse  Hessian  of  the  objective  function. 


32f(P)    32f(P) 


[■ 


V2f  (P) 


1-1 


3i|/' 


3X3^ 


det[v2f(P)]    1     32f(P)  32f  (P) 


3A3i/> 


3X' 


where    the   determinant   of    the    Hessian 


deJv2f(P)l  =  5  f si   risi 


32f  (P) 

3X3^ 


Substituting  the  inverse  Hessian  in  equation  (3-5)  and  expand- 
ing the  equation  by  matrix  multiplication  of  the  inverse 
Hessian  times  the  gradient  yields  a  new  test  point  with  6  the 
only  unknown. 

25 


P.^.    =   P.  -  6 


82f(Pj)     3f(Pj)  92f(Pj)     9f(Pj) 


9ip" 


9X 


dXdip  d\\) 


92f(Pj)     92f(Pj)        /92f(Pj)\2 


3A' 


3ip' 


3X3i|j 


92f(P.)     3f(Pj)  92f(Pj)     9f(Pj) 


3X' 


9tp 


9X9^  9X 


32f(Pj)     92f(Pj)        /92f(Pj)\2 


3X' 


3i/T 


3X3ip 


(3-6) 


The    solution   of   equation    (3-6)    in    turn    requires    the    cal- 
culation  of    all    first    and    second   partial   derivatives    of   the 
objective    function. 


I    n         2* 

3     E     wr 


3f(P)  \i=l      X 


3X 


3X 


n 

E       2p.  (p.B.     +    q. 
1=1 


3B. 


(3-7) 


Similarly, 


9f  (P) 

9ip 


n  9B. 

2       2p.  (p.B.    +   q.    -    6  .)    — — 


(3-8) 


32f  (P) 

3X2 


92f  (P) 


3ijr 


92f  (P) 

3X3i|; 


n 
=      Z 
i=l 


n 
=      E 
i=l 


n 

Z 

i=l 


3B 


2      n 
3*B 


2?i  |»3rl    +  2?i  (PiBi  +  qi-ei)  ^ 


_J 

(3-9) 


/3B. 


2p: 


3^B 


I 


9ip 


-         +    2p.  (p.B.  +  q.  -  0  .  )     ~ 

l       l    l         l         l      «  ,  2 


3\J>' 


3B.     9B 


(3-10) 
32B. 


2P-       VT^   YTi+  2p.  (p.B.   +  q.  -  6  .)     rr^T 
*i      3X       3i|;  ^i^i    l      Mi         l      9X3ip 

(3-11) 
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Solving  equations  (3-7)  through  (3-11)  requires  the 
lengthy  procedure  of  taking  all  first  and  second  partial 
derivatives  of  equation  (2-16),  the  expression  for  angle  B. 
Solutions  to  these  equations  are  contained  in  the  FORTRAN 
Fix  Program  included  in  this  paper. 

C.   GOLDEN  SECTION  SEARCH 

From  equations  (2-15)  and  (3-1) ,  the  objective  function 
can  be  written, 


n  2 

Minimize    E   (p.B.  +  q.  -  8.)  , 
.  ,  L l  l    ^i     l 
1=1 


(3-12) 


the  sum  of  the  squared  azimuth  differences. 

Letting  x  denote  the  indicated  move  in  longitude  and  y 
denote  the  indicated  move  in  latitude  given  by  equation  (3-6)  , 

R  =  (AX,  ^)t 
t 


=  (x,y) 


Th 


en  angle  B.,  given  in  equation  (2-16),  becomes  at  P .  , 


B.  =  cos 

l 


-1 


sin  (^  .  -  Sy)  cos^.  -sin<K  cos  (^  .  -  Sy) 


cos (X .-6x)  cos  A .  +  sin  (X  .  -  5x)  sinX. 
D  i        3  i 


•  r   -1 

in  cos 

n(^.-6y)  sin^.  +  cosCf*.  -  6y)  cos^.  [cos  (X  .  -  6x) 

cosA.  +  sin  (A  .  -  <5x)  sinA. 
i         D  iJ 


(3-13) 


where  subscript  j  denotes  test  point  coordinates  and  sub- 
script i  denotes  DF  site  coordinates.   Since  test  point 
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coordinates  are  arbitrarily  known  and  DF  site  coordinates 
are  constant  for  the  problem,  the  substitution  of  equation 
(3-13)  into  equation  (3-12),  the  objective  function,  results 
in  the  problem  of  minimizing  a  function  with  a  single  unknown, 
6,  the  scalar  coefficient. 

n        2 
Minimize    f(6)    =      Z      oof  (6)  (3-14) 

i=l   1 

The  objective  function  could  be  evaluated  at  incremental 
values  of  6,  say  over  the  interval  (a,b) .   Then  that  6  which 
yields  the  smallest  functional  value  would  be  chosen  as  the 
scalar  coefficient  to  adjust  the  move  indicated  by  the 
second-order  gradient  search.   The  Golden  Section  Search 
(reference  1)  provides  a  fast  and  methodical  approach  to  find 
this  minimizing  6.   Figure  6  illustrates  the  problem;  however, 
the  actual  curve  representing  f (6)  is  unknown  at  the  outset 
of  the  problem. 


\ 

\ 
\ 


\ 


f  (6) 


\  /^f(60) 


\      S 


2 


-o 


6 


«j  S2  b 

Reducing  the  Search  Interval 
FIGURE  6 
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Based  on  the  assumption  that  f (6)  has  a  local  minimum  on 
(a,b) ,  the  objective  is  to  find  that  minimum  without  evalu- 
ating the  function  at  too  many  points  in  (a,b).   The  Golden 
Section  Search  systematically  eliminates  portions  of  the 
interval  by  comparing  functional  values  of  pairs  of  points 
in  (a,b).   From  Figure  6,  f(<5_)  is  less  than  f(6,);  there- 
fore, the  portion  (a, 6,)  can  be  eliminated  in  searching  for 
the  minimum,  reducing  the  search  interval  to  (6,,b).   A  new 
point  then  is  selected  for  comparison  with  f  (6^)  ,  resulting 
in  a  second  elimination  of  whichever  end  interval  corresponds 
to  the  higher  functional  value. 


Letting  I   represent  the  length  of  the  n    search  inter- 


val , 


I.  =  b  -  a 


Two  different  values  of  8    are  tested  at  each  trial.   By 
letting  the  two  values  be  symmetric  on  I  ,  the  length  of 
I   ,  is  the  same  regardless  of  which  end  interval  is  dropped. 
Letting  a  denote  the  relative  position  of  6,  in  L, 


a   =  (61  -  a)/(b  -  a) 


Then 


61  =  a  (b  -  a)  +  a  (3-15) 

Since  6~  is  symmetric  with  6,,  the  relative  position  of  S«  in 
I,  is  given  by  (1  -  a) ,  so  that 

1  -  a  =  (62  -  a)/(b  -  a)  , 
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and 

62  =  (1  -  a)  (b  -  a)  +  a  (3-16) 

Provided  that  6_  yielded  the  minimum  functional  value,  as  in 
Figure  6,  then  6„  becomes  the  first  element  of  the  symmetric 
pair  in  I„  and  its  relative  position  in  I„  is  given  by 

a  =  (<52  -  61)/(1  -  a)  (b  -  a) 

=  \{1   -  o)   (b  -  a)  +  a  -  611/(1  -  a)  (b  -  a) 

from  equation  (3-16) ,  so  that 

61  =  (1  -  a)2(b  -  a)  +  a  (3-17) 

From  equations  (3-15)  and  (3-17) ,  it  follows  that 

a  (b  -  a)  +  a  =  (1  -  a)2(b  -  a)  +  a 

cr  =  (1  -  a)2  (3-18) 

If,  on  the  other  hand,  6,  yields  the  minimum  functional  value, 
then  6,  becomes  the  second  element  of  the  symmetric  test  pair 
in  I?,  and  its  relative  position  in  1^  is  given  by 

1  -  a  =  (61'  -  a)/(l  -  a)  (b  -  a) 

So  that 

6,  =  (1  -  a)2 (b  -  a)  +  a, 

as  in  equation  (3-17)  ,  and  leading  to  the  same  result  as  in 
equation  (3-18) . 

From  equation  (3-18) 

a2  -  3a  +  1  =  0 

which  can  be  solved  by  the  quadratic  formula,  so  that 


30 


a  =    [3  ±  (9  -  4)^1/2 
=  0.382 
Substituting  this  result  in  equations  (3-15)  and  (3-16), 

61  =  0.382  (b  -  a)  +  a  (3-19) 

62  =  0.618  (b  -  a)  +  a  (3-20) 

Generalizing,  let  b   and  a   represent  the  endpoints  of  the 
3       n      n    c  c 

n    search  interval.   Then 

6,  =  0.382  (b   -  a  )  +  a  (3-21) 

1  n    n     n 

6n  =  0.618  (b   -  a  )  +  a  (3-22) 

2  n     n     n 

After  the  first  comparison  on  I, ,    either  6,  or  6  ?  is  the 
survivor  from  the  preceding  iteration,  so  that  a  comparison 
can  be  made  after  evaluating  only  the  new  test  point  symme- 
tric with  the  survivor.   The  Golden  Section  Search  terminates 
when  b   -  a   <  e,  where  e   may  be  tailored  to  the  magnitude  of 
the  indicated  move. 

The  value  of  6  which  minimizes  f(6)  over  (a, ,b, )  is 
entered  in  equation  (3-6)  ,  resulting  in  a  new  test  point 

Pj  +  1  =  (X..,  4'j)t  +  6(x,y)t  (3-23) 


=  P  .  +  6R 


D.   FORTRAN  FIX  PROGRAM 

The  equations  presented  in  this  paper  were  coded  in 
FORTRAN.   The  resulting  program  has  been  successfully 
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executed  on  an  IBM  360/67  computer  system  in  the  W.  R.  Church 
Computer  Center  at  the  Naval  Postgraduate  School. 
1.   Program  Description 

The  program  contains  an  identification  number,  longi- 
tude, latitude  and  name  for  each  DF  site  in  an  imaginary 
network.   Inputs  to  the  program  are  DF  site  identification 
number  and  signal  azimuth.   After  DF  site  names  and  coordi- 
nates are  correlated  with  input  site  identification  numbers, 
necessary  trigonometric  functions  of  site  coordinates  are 
computed  and  stored  in  an  array.   Searching  for  the  minimum 
sum  of  azimuth  differences  begins  from  (0,0)  with  a  Golden 
Section  Search  on  longitude  using  (360,0)  as  the  indicated 
move.   The  result  is  a  minimizing  longitude  on  the  equator 
from  which  a  Golden  Section  Search  is  conducted  on  geocentric 
latitude  with  indicated  move  (0,180).   From  this  minimizing 
geocentric  latitude,  another  Golden  Section  Search  is  con- 
ducted on  longitude  with  (180,0)  as  the  indicated  move. 
Such  alternating  searches  on  longitude  and  geocentric  lati- 
tude continue  until  diminishing  returns  are  indicated  by  a 
new  position  being  less  than  500  nautical  miles  from  the 
previous  minimizing  position.   The  initial  interval  for  the 
Golden  Section  Search  is  set  at  (-1,1),  the  endpoints  for 
equations  (3-19)  and  (3-20).   Subsequent  endpoints  are 
modified  by  the  search  results ,  producing  values  for  equa- 
tions (3-21)  and  (3-22) . 

The  program  then  executes  second-order  gradient 
searches  refined  by  Golden  Section  Searches  on  indicated 
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moves  until  the  most  recent  move  is  less  than  some  input 
termination  value,  say  one-tenth  of  a  nautical  mile.   Indi- 
cated moves,  given  by  equation  (3-6)  ,  are  unbounded  and 
might  occasionally  go  off  the  coordinate  system  of  Figure  1. 
Scaling  any  such  wild  moves  to  the  coordinate  system  by 
modulo  arithmetic  destroys  both  the  direction  and  magnitude 
of  the  vector.   The  indicated  moves  are  kept  in  bounds  by 
conversion  to  a  unit  vector  in  radians  until  such  moves  are 
of  magnitude  less  than  30  nautical  miles,  such  a  small  step 
indicating  that  the  solution  is  at  hand.   By  converting  the 
indicated  move  to  a  unit  vector  prior  to  refinement  by  a 
scalar  coefficient  from  the  Golden  Section  Search,  moves 
from  the  second-order  gradient  search  are  limited  to  about 
3440  nautical  miles  per  iteration. 

Search  output  on  successive  iterations  is  optional, 
depending  on  an  operator  entry.   Unless  search  progress  is 
requested,  searching  terminates  after  thirty  iterations  even 
if  the  latest  step  is  greater  than  the  input  termination 
value.   After  every  tenth  second-order  gradient  search,  if 
such  a  high  number  should  be  required,  the  search  effort  is 
jolted  by  conducting  alternate  Golden  Section  Searches  on 
longitude  with  indicated  move  (90,0)  and  latitude  with  indi- 
cated move  (0,90).   Such  Golden  Section  Search  interruptions 
only  yield  a  new  test  point  when  they  reduce  the  current 
value  of  the  objective  function. 

Prior  to  solution  output,  geocentric  latitude  is 
converted  to  geographical  latitude  using  equation  (2-6) ,  and 


33 


the  solution  is  given  in  degrees  and  minutes  of  longitude, 
east  or  west,  and  latitude,  north  or  south. 
2 .   Results 

The  program  has  worked  successfully  in  over  ten 
thousand  runs  against  simulated  targets  spread  randomly  over 
the  globe  with  various  networks  of  DF  sites  and  various 
numbers  of  sites  intercepting  a  given  target.   In  actual 
practice,  prior  information  about  a  target's  area  of  opera- 
tion might  produce  a  better  starting  point  than  (0,0);  how- 
ever, the  time  required  to  type  such  starting  coordinates 
on  a  keyboard  for  input  would  exceed  the  total  time  required 
for  the  solution  starting  from  (0,0),  which  averages  less 
than  a  second  and  includes  some  diagnostic  output  useful  in 
analyzing  results.   Figure  7  illustrates  a  representative 
set  of  search  iterations  with  Golden  Section  Searches 
plotted  as  dashed  lines  and  second-order  gradient  searches 
solid. 


(0,0) 


Searching 
FIGURE  7 
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IV.   TARGET  LOCATION  CONFIDENCE  PEGION 

That  location  yielding  the  minimum  sum  of  squared  azimuth 
differences,  or  least  squares  solution,  would  coincide  with 
the  location  of  the  signal  emitter  if  all  measured  signal 
azimuths  with  respect  to  the  DF  sites  were  free  from  error. 
Since  then  all  great  circle  arcs  corresponding  to  signal 
azimuths  would  intersect  at  a  point,  the  minimum  value  of 
the  objective  function  would  be  equal  to  zero,  as  is  always 
the  case  when  only  two  signal  azimuths  are  available  to 
compute  a  fix.   Because  the  sites  in  a  DF  network  do  not 
produce  errorless  measurements,  solutions  rarely  coincide 
with  the  emitter  location,  and  there  is  the  question  of  how 
close  the  emitter  might  be  to  the  least  squares  solution. 

A.   SIGNAL  AZIMUTH  ERROR  DISTRIBUTIONS 

Supposing  that  signal  azimuth  measurements  made  at  indi- 
vidual DF  sites  are  not  biased  in  any  of  the  possible  direc- 
tions, then  it  follows  that  signed  errors  have  a  mean  value 
of  zero.   If  such  measurements  were  found  to  be  biased  in 
certain  directions  during  routine  calibration  trials,  then 
the  amount  of  bias  could  be  applied  as  a  correction  to 
subsequent  measurements ,  producing  an  estimated  mean  value 
of  zero.   Lacking  any  reason  to  assume  a  different  distribu- 
tion and  following  the  lead  given  in  the  literature  of 
direction  finding,  the  error  distribution  can  be  taken  to 
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2 

be  normally  distributed  with  mean  zero  and  variance  a., 

assumed  known  and  independent  from  site  to  site. 

The  error  distribution  can  be  made  to  be  symmetric 
around  zero  by  following  the  sign  convertion  that  clockwise 
errors  are  negative  and  counter-clockwise  errors  are  posi- 
tive in  sign. 

B.   THE  OBJECTIVE  FUNCTION  VALUE  AS  A  CHI-SQUARE  STATISTIC 

From  standard  textbooks  on  probability  theory  (for 
example,  reference  2),  it  can  be  verified  that  if  a  random 

variable  is  normally  distributed  with  mean  zero  and  variance 

2 

a  ,  then  the  square  of  that  random  variable  is  governed  by 

the  Gamma  distribution  with  parameters  (l/2,l/2a  ).   Further- 
more, the  sum  of  n  independent  Gamma  distributed  random 
variables  divided  by  the  variance  of  the  underlying  normal 
distributions  is  a  Chi-Square  statistic  with  n  degrees  of 
freedom.   Hence,  where  co  denotes  the  azimuth  difference, 

Uj^NtOrO?)  (4-1) 

.f   J  "N^ln)  (4"3) 

i 

Since  the  normal  distribution  here  referred  to  applies 
to  the  azimuth  differences  between  the  measured  signal 
azimuth  and  the  azimuth  of  the  unknown  emitter  location,  the 
co.  refer  to  actual  and  unknown  errors  rather  than  calculated 
differences  between  solution  azimuths  and  signal  azimuths. 
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Equation  (4-3)  therefore  specifies  the  distribution  of 
values  of  the  objective  function  (equation  (3-12) )  at  unknown 
emitter  locations  under  the  assumption  that  azimuth  errors 
are  normally  distributed  with  known  variance. 

C.   CONFIDENCE  LEVELS  AND  REGION  BOUNDARIES 

Since  the  sum  of  squared  azimuth  errors  divided  by 
variances  is  a  Chi-Square  statistic  under  the  previous 
assumptions,  there  is  a  specific  probability/  denoted  a, 
that  this  sum  will  be  less  than  or  equal  to  any  given  value. 
Equivalently ,  for  a  fixed  probability  a,  or  confidence 
level,  there  corresponds  an  upper  bound,  denoted  C,    *  , 
which  depends  upon  the  degrees  of  freedom  as  well  as  the 
confidence  level. 

I        4r  OK   <  C,     J 

i=l  o\      x  ~   (n'a) 

The  choice  of  confidence  level,  a,  is  arbitrary,   The 

corresponding  value  for  C  ,    s  can  be  found  in  Chi-Square 
c  3  (n,a)  1 

tables . 

For  a  given  fix  problem,  the  value  of  the  sum  of  least 
squares,  denoted  f*,  is  the  minimum  possible  value  of  the 
Chi-Square  statistic  for  the  observed  signal  azimuths. 
Furthermore,  by  evaluating  the  objective  function  at  suc- 
cessive points  outward  in  any  direction  from  the  least 
squares  solution,  it  can  be  seen  that  the  farther  away  the 
point  evaluated  up  to  the  antipode ,  the  higher  the  objective 
function  value.   The  location  corresponding  to  f*  can 
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therefore  be  regarded  as  a  point  in  a  region  which  contains 

all  the  possible  emitter  locations  for  which  the  value  of 

the  objective  function  is  less  than  or  equal  to  an  upper 

bound  C,    ,  ,  given  the  set  of  observed  signal  azimuths, 
(n  ,a)   ^ 

Such  a  region  is  bounded  by  all  locations  for  which  the 

objective  function  has  a  value  equal  to  C,    , ,  and  it  can 

^         (n,a) 

be  stated  that,  for  the  observed  signal  azimuths,  there  is  a 
probability  a  that  the  bounded  region  will  contain  the 
unknown  emitter  location. 

Daniels  and  Vajda  (1951)  challenge  the  utility  of  such 
an  approach  on  the  grounds  that  at  given  confidence  levels 
the  corresponding  region  may  not  exist  with  positive 
probability  if  n  is  greater  than  two.   Their  observation  is 
a  consequence  of  the  fact  that  except  for  n  equal  to  two 
all  of  the  great  circle  arcs  corresponding  to  signal  azimuths 
seldom  intersect  at  a  single  point   thereby  yielding  f* 
equal  to  zero. 

Prior  to  the  availability  of  the  high-speed  digital 
computer,  the  practical  problems  of  producing  such  a  region 
were  certainly  more  striking  than  theoretical  objections. 
Even  so,  the  Daniels/Vajda  criticism  was  more  than  casual, 
which  may  indicate  a  requirement  for  additional  arguments 
to  justify  the  use  of  such  a  confidence  region. 

Suppose  that  f*  is  equal  to  k  for  some  least  squares 
solution,  where 


(4-5) 


and  the  probability  p  may  be  determined  by  interpolating 

with  the  value  k  in  a  table  of  Chi-Squared  values.   Daniels 

and  Vajda  are  concerned  that  if  the  confidence  level  a  had 

been  initially  specified  at  some  value  less  than  p,  then 

there  exists  no  region  of  locations  having  objective  function 

values  less  than  or  equal  to  k ,  k  being  the  minimum  possible 

value  on  the  globe  for  the  given  signal  azimuth  observations. 

There  are  at  least  two  responses  to  such  an  occurrence  and 

the  accompanying  empty  region  objection:   (1)  Increase  the 

confidence  level  whenever  the  region  does  not  exist  at  the 

prior  confidence  level  so  that  C,    s  is  greater  than  k,  and 
r  (n,a)     r 

(2)  Routinely  specify  a  confidence  level  of  sufficient  mag- 
nitude to  make  the  occurrence  of  an  empty  region  highly 
improbable. 

Either  of  the  above  alternatives  insures  the  production 
of  a  non-empty  confidence  region.   There  remains  the  question 
of  what  is  represented  by  a  possible  contention  inferable 
from  equation  (4-5) ,  namely  that  positive  probability  p  is 
attached  to  a  single  point,  the  least  squares  solution.   Any 
such  statement  would  be  puzzling  when  approached  with  the 
viewpoint  that  specific  random  variables  on  a  continuum  have 
zero  probability  and  that  probability  statements  can  only  be 
made  about  boundaries  for  continuous  random  variables  such 
as  location  coordinates. 

An  answer  may  be  attempted  by  first  pointing  out  that 
azimuth  errors  are  unknown  and  independent  random  variables, 
but  they  are  attached  to  actual  and  unknown  emitter  azimuths 
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with  respect  to  the  set  of  intercepting  DF  sites.   Observed 
signal  azimuths  have  two  components,  the  first  of  which, 
emitter  azimuth,  is  dependent  upon  emitter  location.   The 
second  component,  azimuth  error,  is  independent  from  site  to 
site  and  is  the  theoretical  basis  for  the  Chi-Square  contour, 
confidence  region.   Prior  to  the  production  of  a  particular 
confidence  region,  the  normality  assumption  permits  the 
specification  of  upper  bounds  for  the  Chi-Square  statistic 
corresponding  to  various  levels  of  confidence.   During 
production  of  the  particular  confidence  region,  the  require- 
ment that  a  single  location  must  satisfy  the  observations 
places  a  lower  bound  on  the  Chi-Square  statistic. 

Each  emitter  azimuth  can  depart  from  the  observed  signal 
azimuth,  but  there  are  conceptual  bands  (based  .on  the  error 
distribution  assumptions)  bounding  such  departures  for  any 
given  probability.   When  more  than  two  such  bands  are  over- 
lapped in  the  solution  region,  then  a  proportion  of  the 
possible  combinations  of  such  departures  are  ruled  out  by 
the  constraint  that  azimuth  errors  are  attached  to  emitter 
azimuths.   Whenever  f*  is  greater  than  zero,  then  the  prob- 
ability, denoted  p* ,  associated  with  f*  as  an  upper  boundary 
from  equation  (4-5)  represents  the  proportion  of  impossible 
sums  of  squares  out  of  those  permitted  by  the  density  func- 
tions governing  azimuth  errors. 

A  thought  experiment  may  clarify  these  remarks.   Suppose 
it  were  possible  to  specify  all  combinations  of  azimuth 
errors  permitted  by  their  respective  density  functions  prior 
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to  finding  the  least  squares  solution.   There  exists  a  loca- 
tion satisfying  one  of  the  postulated  combinations  of  azimuth 
errors  for  a  specific  set  of  signal  azimuths  only  if  applica- 
tion of  each  of  the  postulated  azimuth  errors  in  the  combina- 
tion to  its  corresponding  signal  azimuth  yields  a  corrected 
set  of  emitter  azimuths  whose  corresponding  great  circle  arcs 
intersect  at  a  point.   Few  such  postulated  combinations  in 
the  infinite  set  can  be  satisfied  since  the  point  of  inter- 
section for  any  two  corrected  emitter  azimuths  uniquely 
determines  the  required  azimuth  corrections  for  the  rest  of 
the  observations.   Furthermore,  there  cannot  exist  a  satis- 
fying location  for  any  postulated  combination  of  azimuth 
errors  having  a  sum  of  squares  less  than  f*  since  no  such 
corrections  are  of  sufficient  magnitude  to  result  in  a 
common  intersection. 

These  suggestions  can  be  illustrated  by  an  extreme  and 
improbable  example  of  emptiness  of  the  confidence  region  for 
three  signal  azimuth  great  circle  arcs  forming  a  cocked  hat. 
Suppose  that  the  relationships  between  the  DF  sites  and  the 
emitter  are  as  given  in  Figure  8,  and  suppose  that  the  errors 
associated  with  the  three  DF  sites  are  all  governed  by  the 
standard  normal  distribution  (variance  equal  to  one  degree 
squared) .   Suppose  also  that  each  of  the  three  signal 
azimuths  were  in  error  by  minus  two  degrees  (where  the 
probability  that  azimuth  error  is  less  than  or  equal  to 
minus  two  degrees  is  equal  to  .0227)  in  order  to  form  one 
of  the  two  largest  possible  cocked  hats  for  [w. |  equal  to 
two  degrees.   Then  the  least  squares  solution  will  be  found 
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near  the  center  of  the  cocked  hat,  and  f*  is  about  equal  to 
twelve  degrees  squared.   For  three  degrees  of  freedom,  there 
was  a  .992  probability  that  the  Chi-Square  statistic  would 
be  less  than  or  equal  to  twelve  degrees  squared,  indicating 
an  empty  confidence  region  at  the  .992  level  (or  any  lower 
levels) .   Now  the  probability  that  |w. |  is  greater  than  four 
degrees  is  equal  to  .000,  suggesting  that  a  boundary  of  four 
degrees  on  either  side  of  each  signal  azimuth  great  circle 
arc  should  contain  all  possible  emitter  azimuths  from  which 
each  signal  azimuth  represents  a  departure.   The  situation  is 
illustrated  in  Figure  8. 

With  each  emitter  azimuth  free  to  be  anywhere  within  its 
dashed  band  by  the  normality  assumption,  but  dependent  upon 
the  actual  emitter  location  and  therefore  constrained  by  the 
other  dashed  bands ,  the  shaded  intersection  of  the  three 
bands  must  (for  any  non-negligible  probability)  contain  the 
emitter.   Moreover,  the  Chi-Square  contour  methodology  would 
produce  a  circular  confidence  region  within  the  shaded 
triangle  for  any  confidence  level  greater  than  .992.   The 
center  of  the  shaded  region  corresponds  to  the  least  squares 
solution.   Suppose  a  "center"  to  be  the  least  squares  solu- 
tion united  with  all  points  forming  a  neighborhood  within 
radius  e    of  the  center.   Then  this  "center"  had  a  probability 
of  .992  +  f(e)  of  containing  the  emitter  location.   Thus 
there  was  a  probability  of  .008  -  f (z)    of  finding  the  emitter 
within  a  region  bounded  by  a  "centerless"  circle  inscribed 
in  the  shaded  triangle.   There  is  zero  probability  that  the 
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emitter  lies  outside  such  an  inscribed  circle,  or  confidence 
region  at  the  100  percent  level.   The  description  of  the 
confidence  region  as  circular  applies  to  Figure  8  only 
because  the  DF  sites  are  equidistant  and  symmetric  around 
the  emitter,  and  because  azimuth  errors  are  minus  two  degrees, 

The  implication  of  the  foregoing  is  that  if  errors  are 
normally  distributed  with  known  variance,  then  abnormally 
large  cocked  hats  would  be  indicative  of  excellent  solutions 
via  the  least  squares  approach.   The  point  would  be  difficult 
for  operators  to  accept  since  high  confidence  is  normally 
associated  with  a  set  of  great  circle  arcs  corresponding  to 
signal  azimuths  that  nearly  intersect  at  a  point,  and  large 
departures  are  suspect.   The  discrepancy  between  the  above 
conclusion  and  experience  may  be  attributed  to  the  fact  that 
abnormally  large  cocked  hats  only  occur  when  outliers  (signal 
azimuth  errors  that  do  not  obey  the  empirically  estimated 
distribution)  enter  the  data  base. 

Sometimes  strange  relationships  between  the  sizes  of 
cocked  hats  and  accuracy  of  locations  (least  squares  solu- 
tions) can  also  be  demonstrated  with  a  network  of  three  DF 
sites  with  errors  governed  by  a  Bernoulli  distribution, 
having  equal  probability  that  azimuth  error  is  plus  or  minus 
three  degrees.   Figure  9  illustrates  the  situation.   Given 
an  emitter  at  the  center  of  the  star,  observed  signal  azi- 
muths are  in  error  by  plus  or  minus  -three  degrees ,  each 
with  probability  1/2.   All  possible  configurations  of 
cocked  hats  are  therefore  represented  on  the  figure  for  the 
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sketched  relationship  between  DF  sites  and  emitter  location. 
Six  of  the  possible  cocked  hats  are  the  small  triangles  at 
the  tips  of  the  star.   The  other  two  possibilities  are  the 
union  of  the  body  of  the  star  with  the  small  triangles  at 
points  NE ,  SE  and  W,  and  the  union  of  the  body  of  the  star 
with  the  small  triangles  at  points  E,  SW  and  NW.   The  latter 
two  cocked  hats  are  abnormally  large  and  have  f*  equal  to  27 
degrees  squared  (which  is  over  twice  the  value  of  a  Chi-Square 
statistic  at  the  .995  level  for  a  normal  error  distribution). 
Yet,  the  center  of  these  largest  cocked  hats  coincides 
approximately  with  the  emitter  location,  indicating  excellent 
least  squares  solutions ,  and  the  centers  of  the  triangles  at 
the  points  of  the  star,  the  smaller  cocked  hats,  are  far 
from  the  emitter  location,  even  though  their  f*  value  is 
relatively  small. 

The  probability  of  occurrence  of  such  large  cocked  hats 
when  error  is  governed  by  a  normal  rather  than  Bernoulli 
distribution  is  negligible,  and  they  do  not  appear  when  data 
is  simulated  with  a  normal  error  distribution  having  reason- 
able variance.   However,  in  about  400  simulations  during 
development  of  the  Fix  Program  they  were  forced  to  appear 
by  governing  error  with  such  a  Bernoulli  distribution.   As 
would  be  expected  from  Figure  9,  a  plot  of  location  error 
versus  f*  verified  that  location  error  decreases  as  f* 
approaches  its  maximum  of  27  degrees  squared  for  this  unusual 
error  distribution.   Such  a  simulation  naturally  yields  about 
25  percent  of  the  abnormally  large  cocked  hats  and  75  per- 
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cent  triangles  at  the  points  of  the  star.   With  a  random 
distribution  of  DF  sites  and  signal  emitters  over  the  globe, 
the  stars  are  not  so  well-shaped  as  that  in  Figure  9. 

If  one  regards  the  fix  problem  with  least  squares  solu- 
tion and  confidence  region  construction  at  confidence  level 
a  as  a  Bernoulli  trial  with  success  if  the  confidence  region 
contains  the  emitter  location  and  failure  if  the  emitter  is 
outside  the  confidence  region,  then  the  expected  number  of 
successes  for  m  independent  trials  is  am.   That  the  Chi- 
Square  contour  confidence  regions  have  normalized  frequency 
of  containment  about  equal  to  the  probability  of  containment 
is  therefore  expected  to  be  true  and  was  demonstrated  to  be 
true  by  simulation.   At  confidence  level  a,  on  about  a  per- 
cent of  the  trials  non-empty  confidence  regions  are  produced 
which  contain  the  actual  emitter.   On  the  remaining  (1  -  a) 
percent  of  the  trials,  outcomes  are  distributed  among  empty 
regions  and  non-empty  regions  with  the  actual  emitter  out- 
side the  confidence  region.   The  lower  the  confidence  level 
the  higher  the  proportion  of  empty  regions  among  the  failures. 

D.   FORTRAN  CONFIDENCE  REGION  PROGRAM 

A  plot  subroutine  is  included  in  the  Fix  Program  of  this 
paper.   Sample  plots  are  included  in  the  Computer  Output 
section  of  the  paper.   The  plot  is  a  square  grid  of  constant 
dimension,  scaled  in  nautical  miles  to  contain  the  confidence 
region.   An  asterisk  on  the  grid  indicates  that  the  corres- 
ponding point  has  an  objective  function  value  greater  than 
the  Chi-Square  boundary  value  for  the  specified  confidence 
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level.   The  least  squares  solution  is  at  the  center  of  the 
grid  and  is  indicated  by  the  letter  S.   During  program 
development,  it  was  found  that  a  90  percent  confidence  level 
is  sufficient  to  preclude  the  occurrence  of  empty  confidence 
regions,  and  the  sample  regions  are  of  that  size. 
1.   Program  Description 

Given  in  a  data  statement  of  the  subroutine  are 
Chi-Square  values  corresponding  to  the  numbers  of  possible 
degrees  of  freedom,  up  to  n  for  a  network  of  n  DF  sites,  and 
the  selected  confidence  level.   The  confidence  region  boun- 
dary value  is  selected  according  to  the  number  of  signal 
azimuths  available  to  find  a  solution. 

The  dimension  of  the  grid  is  41  rows  and  51  columns, 
tailored  to  the  line  printer  distances  between  rows  and 
columns,  so  that  the  ratio  of  the  distances  between  rows  and 
columns  is  4:5.   Depending  upon  the  scale  of  the  grid,  coor- 
dinates on  the  globe  are  calculated  to  match  each  of  the 
2091  points  on  the  grid.   Curvature  of  the  earth  calculations 
are  used  to  determine  distances  between  points  on  the  grid, 
thereby  avoiding  distortion  inherent  in  planar  assumptions. 

At  each  of  the  points  on  a  row,  the  sum  of  squared 
azimuth  differences  is  computed  for  comparison  with  C,    ,,' 
an  asterisk  being  assigned  for  grid  entry  if  C,    x  is 
exceeded.   After  the  confidence  region  is  exited  on  a  given 
row,  remaining  points  on  the  row  are  assigned  asterisks 
without  computation.   Similarly,  after  the  confidence  region 
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is  completed,  the  rest  of  the  grid  is  filled  with  asterisks 

without  computation. 

The  area  of  the  confidence  region  is  estimated  by 

counting  the  blanks  and  multiplying  their  sum  by  the  area 

represented  by  each  cell,  80  square  miles  in  the  case  of 

a  400  by  400  nautical  mile  grid.   If  computations  on  the 

first  row  yield  any  points  with  objective  function  value  less 

than  or  equal  to  C,    N ,  then  the  grid  is  rescaled  to  cover 
^         (n,a)  * 

a  greater  area,  and  the  computation  starts  over.   If  the 
confidence  region  leaks  out  on  either  of  the  sides  of  the 
grid  during  computation,  then  again  the  grid  is  rescaled, 
and  computation  starts  over.   In  practice,  the  initial  area 
covered  by  the  grid  can  be  tailored  to  the  network  and 
expected  size  of  confidence  regions. 
2 .   Results 

If  the  solution  search  in  the  Fix  Program  should 
terminate  in  the  region  of  the  least  squares  solution  before 
the  least  squares  point  is  found,  then  the  confidence  region 
is  noticeably  asymmetric  about  the  center  at  S.   The  lack  of 
symmetry  indicates  that  the  fix  solution  should  be  recomputed 
with  a  smaller  termination  step.   The  center  of  the  confi- 
dence region,  or  least  squares  location,  could  also  be 
estimated  from  the  grid  itself,  but  only  to  an  accuracy 
limited  by  the  dimensions  of  a  cell. 

The  confidence  region  subroutine  was  used  successfully 
in  over  a  thousand  runs  with  simulated  data.   Its  chief 
deficiency  is  the  computer  time  required  to  produce  the  grid 
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entries,  averaging  slightly  over  ten  seconds  for  numbers  of 
signal  azimuths  per  target  varying  uniformly  from  three  to 
nine.   On  a  slower  machine  than  the  IBM  360/6  7,  the  subrou- 
tine would  require  revision  to  be  usable  for  production. 
Sampling  every  other  or  every  third  cell  would  be  possible 
avenues  for  reducing  computation  time,  but  the  present 
version  was  preferred  for  a  more  accurate  estimate  of  con- 
fidence region  area.   A  more  sophisticated  method  for  reduc- 
ing computation  time  would  be  to  use  a  modified  Golden 
Section  Search  to  find  the  confidence  region  boundary  outward 
from  the  center  of  the  region,  store  coordinate  pairs  of  the 
boundary  in  an  m  by  n  array,  and  then  print  the  boundary  on 
a  scaled  grid. 
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V.   INTERCEPT  GEOMETRY 

The  Fix  Program  produces  several  statistics  that  are 
useful  in  describing  intercept  geometry.   The  first  group 
of  these  statistics  is  related  to  signal  azimuths  and  is 
used  to  plot  position  lines  as  linear  approximations  of  the 
great  circle  arcs  in  the  region  of  the  solution.   In  the 
second  group  are  some  values  to  which  location  error  is 
sensitive. 

A.   PLOTTING  POSITION  LINES  IN  SOLUTION  REGION 

Following  the  optional  summary  of  search  progress,  the 
least  squares  solution  and  the  value  of  f*,  the  computer 
output  has  a  table  of  values  in  which  each  row  corresponds 
to  one  of  the  intercepting  DF  sites.   Row  entries  include 
the  site  identification  number,  the  site  name,  the  azimuth 
of  the  solution,  the  signal  azimuth,  their  difference,  the 
range  in  nautical  miles  from  the  solution  to  the  site,  and 
three  values  used  to  plot  the  position  line  relative  to  the 
least  squares  solution.   These  values  are  calculated  by 
analogy  with  a  right  spherical  triangle  illustrated  in 
Figure  10. 

Figure  10  is  labelled  to  correspond  to  the  base  of  the 
polar  triangle  of  Figure  3  with  the  DF  site  at  vertex  B  and 
the  least  squares  solution  at  Vertex  C.   Vertex  angle  A  is 
a  right  angle  formed  by  connecting  C  to  the  closest  point  of 
approach  (CPA)  of  the  great  circle  arc  related  to  the  signal 
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Right  Spherical  Triangle 
FIGURE  10 

azimuth  6.   Arc  a,  connecting  the  DF  site  and  least  squares 
solution,  is  related  to  the  solution  azimuth,  denoted  3*/ 
and  is  given  by  equation  (2-11) ,  using  the  polar  triangle  of 
Figure  3.   By  the  lav;  of  sines  for  spherical  triangles, 


sin  b  =  sinoo  sin  a 


(5-1) 


The  curvature  of  earth  distance  covered  by  arc  b  is 
closely  approximated  by  multiplying  arc  b  in  radians  by  the 
mean  earth  radius  of  3439.9  nautical  miles,  yielding  the 
nautical  mile  departure  of  the  CPA  from  the  solution. 

The  site  azimuth  with  respect  to  the  solution,  denoted 
y* ,    is  found  after  calculating  angle  C  of  Figure  3  in  the 
same  manner  that  vertex  angle  B  and  azimuth  3  were  calculated 
in  equations  (2-12)  and  (2-13) .   Then  the  azimuth  of  the  CPA 
with  respect  to  the  solution  is  the  sum  of  y*    and  C,  in  the 
case  illustrated.   Finally,  the  azimuth  of  the  DF  site  with 
respect  to  the  CPA  is  found  by  subtracting  90  degrees  from 
the  sum  of  y*and  C,  in  the  case  illustrated.   Similar  argu- 
ments apply  to  three  other  general  cases,  solution  east  or 
west  of  the  DF  site  and  w  positive  or  negative,  and  resulting 
conditional  equations  are  included  in  the  Fix  Program. 
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The  azimuth  and  distance  of  the  CPA  with  respect  to  the 
solution  are  the  components  of  a  vector  which  can  be  plotted 
on  polar  coordinate  paper.   A  line  of  position  perpendicular 
to  this  vector  approximates  the  great  circle  arc  in  the 
region  of  the  solution.   The  direction  of  the  intercepting 
DF  site  is  given  by  the  azimuth  value  in  the  last  column  of 
the  output  row.   The  plot  can  also  be  constructed  with  a 
parallel  ruler  by  orienting  the  ruler  along  the  site  azimuth 
with  respect  to  the  CPA,  moving  the  free  edge  the  distance 
to  scale  of  the  CPA  toward  the  azimuth  of  the  CPA.   Such 
lines  of  position  are  plotted  to  the  grid  scale  on  several 
of  the  confidence  region  plots  provided  in  the  computer  out- 
put section  of  this  paper. 

B.   MEAN  ACUTE  ANGLE  OF  INTERSECTION 

There  are  factors  that  may  be  said  to  reduce  location 
error,  given  that  signal  azimuth  error  is  present  but  un- 
known.  One  such  factor  is  the  number  of  DF  sites  providing 
signal  azimuths  to  make  a  fix:   the  more  independent  measure- 
ments, the  better  the  emitter  location  estimate  by  least 
squares  solution.   Another  factor  is  the  layout  of  the  DF 
network,  the  obvious  supposition  being  that  a  network  spread 
over  a  large  area  produces  better  locations  of  globally 
distributed  targets  than  a  network  confined  to  a  small  area. 

A  slightly  less  obvious  statistic  is  the  magnitude  of 
the  acute  angle  of  intersection  of  great  circle  arcs  corres- 
ponding to  signal  azimuths.   For  two  DF  sites  participating, 
the  dependence  of  location  error  on  the  acute  angle  of 
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intersection  is  illustrated  in  Figure  11.   In  case  a,  two 
lines  of  position  intersect  at  about  90  degrees;  in  case  b, 
the  acute  angle  of  intersection  is  about  ten  degrees.   The 
difference  in  the  sizes  of  regions  contained  by  the  dashed 
lines  demonstrates  the  effect  of  two  degree  uncertainty  of 
measurement  on  either  side  of  a  signal  azimuth  on  the  two 
different  intersections. 

Although  the  confidence  regions  described  in  Section  IV 
of  this  paper  do  convey  this  dependence,  that  for  case  a 
being  a  near  circle  contained  in  the  dashed  lines  and  that 
for  case  b  being  an  elongated  ellipse  of  orientation  and 
dimension  similar  to  to  the  dashed  region,  the  acute  angle 
of  intersection  statistic  was  found  to  be  significant  enough 
to  warrant  isolated  analysis.   It  can  be  demonstrated 
graphically  that  the  dependence  also  exists  for  n  greater 
than  two.   For  all  cases,  the  dependence  can  be  examined  by 
calculating  the  mean  value  of  all  the  acute  angles  of  inter- 
section formed  by  all  possible  permutations  of  two  lines  of 
position,  given  by  (n)?.   The  mean  acute  angle  of  inter- 
section is  part  of  the  Fix  Program  output.   Dependence  of 
location  error  on  this  statistic  is  examined  in  later  sec- 
tions of  this  paper.   For  individual  pairs,  the  acute  angle 
of  intersection  is  closely  approximated  by  taking  differences 
from  the  site  azimuths  with  respect  to  the  CPA. 
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a.   90  Degree  Acute  Angle  of  Intersection 


b.   10  Degree  Acute  Angle  of  Intersection 


FIGURE  11 
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C.   MEAN  ARC  STATISTIC 

A  more  subtle  factor  affecting  location  error  is  a 
function  of  arc  length  or  distance  from  the  DF  site  to  the 
emitter  location.   The  arc  from  the  DF  site  to  the  emitter 
and  all  other  arcs  representing  possible  departures  due  to 
erroneous  signal  azimuth  measurements  converge  at  the  DF 
site  and  at  its  antipode.   Therefore,  the  maximum  nautical 
mile  departure  of  a  line  of  position  from  the  emitter  loca- 
tion for  constant  azimuth  error  occurs  at  90  degrees  of  arc, 
or  about  5403  nautical  miles,  from  the  DF  site. 

The  arc  from  the  DF  site  to  the  least  squares  solution 
is  given  by  equation  (2-11) .   The  arc  statistic  is  equal  to 
this  arc,  or  to  its  supplement  if  the  arc  is  greater  than 
90  degrees.   The  dependence  of  location  error  on  the  arc 
statistic  is  illustrated  in  Figure  12,  where  the  solid  line 
represents  an  arc  from  the  DF  site  through  the  emitter  loca- 
tion to  the  DF  site  antipode,  and  the  dashed  lines  represent 
three  degree  departures  from  this  arc. 

The  distance  from  the  site-to-emitter  arc  to  the  depar- 
ture arcs  increases  non-linearly  from  zero  at  the  DF  site 


Site  Antipode 

Relating  Error  to  Emitter  Distance" 
FIGURE  12 
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and  its  antipode  to  180  nautical  miles  at  the  midway  point. 
In  a  later  section  of  this  paper,  location  error  sensitivity 
to  the  mean  arc  statistic  for  all  DF  sites  participating  in 
the  location  is  examined. 
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VI.   SIMULATION  PROCEDURE 

A  modified  version  of  the  Fix  Program  was  used  to  simu- 
late operation  against  targets.   Targets  were  generated  by 
a  subroutine  which  provided  longitude  and  geocentric  lati- 
tude describing  random  locations  uniformly  distributed  over 
the  globe.   Such  coordinates  were  saved  for  later  comparison 
with  the  least  squares  solution. 

The  number  of  intercepting  sites  for  each  target  was 
either  specified  for  simulations  on  fixed  degrees  of  freedom, 
or  randomly  generated  according  to  a  uniform  distribution 
from  three  to  nine.   In  both  cases,  the  specific  intercept- 
ing DF  sites  were  randomly  selected  from  a  set  of  19  DF  sites 
composing  the  particular  network  being  simulated.   Azimuths 
of  the  target  location  with  respect  to  the  participating  DF 
sites  were  then  calculated.   For  each  such  emitter  azimuth, 
a  random  error  was  added,  with  direction  and  magnitude  of 
the  error  governed  by  the  standard  normal  distribution.   The 
emitter  azimuth  with  error  was  then  treated  as  a  measured 
signal  azimuth  for  input  with  the  site  identification  number 
to  the  Fix  Program. 

During  program  development,  the  simulation  procedure  was 
used  with  zero  signal  azimuth  error  generated  to  insure  that 
the  least  squares  solution  of  the  search  procedure  coincided 
with  the  target  coordinates  to  at  least  an  accuracy  of  three 
decimal  places.   A  search  termination  value  of  one-tenth  of 
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a  nautical  mile  was  found  to  be  adequate,  and  the  present 
combination  of  Golden  Section  Searches  and  second-order 
gradient  searches  was  found  to  be  satisfactory  during  the 
development  simulations.   The  Fix  Program  has  also  been 
successfully  tested  with  up  to  19  participating  DF  sites  on 
a  single  target,  with  a  set  of  DF  sites  contained  in  a  180 
by  240  nautical  mile  area,  and  with  higher  variance  error 
distributions,  including  outliers. 

A.   GEOCENTRIC  AND  GEOGRAPHICAL  LATITUDE 

A  comparison  of  least  squares  solutions  using  geocentric 
and  geographical  latitude  in  Fix  Program  calculations  was 
made  in  order  to  examine  the  accuracy  gain  of  the  geocentric 
model.   Errorless  signal  azimuths  of  random  targets  were 
computed  using  geocentric  latitude  in  calculations.   A  least 
squares  solution  was  then  found  by  using  geographical  lati- 
tude in  calculations,  and  the  conversion  of  the  solution  to 
geographical  latitude  was  omitted.   This  distance  between 
the  geographical  coordinates  of  the  target  and  the  pseudo- 
geographical  coordinates  of  the  least  squares  solution 
measures  the  departure  from  the  least  squares  location  due 
to  use  of  geographical  latitude  on  an  assumed  sphere. 

Radial  departures  for  two  degrees  of  freedom  were  occa- 
sionally devastating  (424  nautical  miles  in  one  case) .   Using 
an  imaginary  network,  called  the  national  network,  and  degrees 
of  freedom  uniformly  distributed  from  three  to  nine,  the  mean 
radial  departure  for  14  8  samples  was  6.8  nautical  miles. 
The  direction  of  such  departures  tends  toward  the  equator  for 
targets  in  both  northern  and  southern  hemispheres. 
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VII.   LOCATION  ERROR  ANALYSIS  WITH  FIXED  DEGREES  OF  FREEDOM 

The  sensitivity  of  location  error  to  several  factors  was 
examined  in  simulations  v;ith  degrees  of  freedom  fixed  from 
two  to  seven.   From  10  0  to  200  samples  were  generated  at  each 
level  for  each  of  the  networks.   All  of  the  following  plots 
are  drawn  to  the  same  scale  to  simplify  between  plot  compari- 
sons of  the  effects  of  different  degrees  of  freedom  and 
different  network  configurations.   Values  too  large  for  the 
common  scale  are  entered  in  a  row  at  the  upper  limit  of  the 
plot  above  an  off  of  scale  indicator. 

A.   COMPARISON  OF  ERROR  ON  NATIONAL  AND  GLOBAL  NETWORKS 

Member  sites  of  the  national  network  varied  in  longitude 
from  Attu  to  Boston  and  in  latitude  from  Anchorage  to  Hono- 
lulu.  The  specific  sites  are  identified  in  the  second  data 
statement  of  the  Fix  Program.   Sites  in  the  widely  distri- 
buted global  network  included  Anchorage,  Bombay,  Boston, 
Buenos  Aires,  Cape  Town,  Chris tchurch,  Dakar,  Honolulu,  Houston, 
Lima,  Naples,  Papeete,  Perth,  Recife,  Seattle,  Seychelles, 
Tehran,  Tokyo  and  Wake. 

Figure  13  is  a  plot  of  mean  location  error  from  each  of 
the  sample  groups  against  the  number  of  intercepting  sites. 
Dashed  lines  connect  the  means  of  each  network.   The  plot 
shows  nearly  double  the  location  error  in  the  national  net- 
work even  at  higher  degrees  of  freedom. 
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With  all  other  factors  random,  the  average  accuracy  of 
a  network  is  dependent  upon  the  network  layout  and  the 
variance  of  signal  azimuth  measurement  error.   For  a  given 
network,  the  mean  location  error  would  vary  in  the  same 
direction  as  changes  in  such  variance,  in  all  of  the  follow- 
ing simulations  held  constant  at  one  degree. 

B.   ERROR  VERSUS  MEAN  ACUTE  ANGLE  OF  INTERSECTION 

Based  on  samples  from  the  national  network,  Figures  14, 
15  and  16  are  scatter  plots  of  location  error  versus  the 
mean  acute  angle  of  intersection  explained  in  section  V-B. 
Degrees  of  freedom  were  fixed  at  three,  five  and  seven, 
respectively. 

The  strong  correlation  in  the  national  network  can  be 
interpreted  as  a  reduction  of  both  the  mean  and  variance  of 
location  error  as  intersections  of  great  circle  arcs  are  less 
acute.   Rough  scatter  plots  for  two,  four  and  six  degrees  of 
freedom  showed  the  same  correlation,  naturally  most  pro- 
nounced at  two  degrees  of  freedom.   The  set  of  plots  shows 
a  shrinking  of  the  range  of  the  mean  acute  angle  of  inter- 
section statistics  as  degrees  of  freedom  are  increased. 

Figures  17,  18  and  19  are  similar  scatter  plots  with 
entries  from  simulations  on  the  global  network.   The  corre- 
lation is  still  evident  at  three  degrees  of  freedom  but  be- 
comes masked  at  higher  levels  as  the  layout  of  the  network 
reduces  the  chance  of  producing  a  combination  of  target  and 
intercepting  sites  with  a  small  mean  acute  angle  of  inter- 
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section.   Consequently,  shrinking  of  the  range  of  the 
statistic  is  extremely  evident  on  Figures  18  and  19. 

C.   ERROR  VERSUS  MEAN  ARC  STATISTIC 

Figures  20,  21  and  22  are  scatter  plots  of  national 
network  location  error  versus  the  mean  arc  statistic 
(explained  in  section  V.C)  for  three,  five  and  seven  degrees 
of  freedom,  respectively.   At  all  levels,  the  large  location 
errors  tend  to  be  associated  with  the  high  end  of  the  range 
of  mean  arc  statistics.   The  strength  of  the  correlation 
appears  to  increase  as  degrees  of  freedom  increase  on 
successive  plots. 

Figures  23,  24  and  25  are  similar  plots  for  the  global 
network.   Again,  the  correlation  is  masked  after  three 
degrees  of  freedom  with  the  shrinkage  effect  extremely  pro- 
nounced as  the  degrees  of  freedom  increase. 
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IIX.   LOCATION  ERROR  ANALYSIS  WITH  VARYING  DEGREES  OF  FREEDOM 

Simulations  were  conducted  on  both  networks  with  degrees 
of  freedom  generated  uniformly  from  three  to  nine.   Chi- 
square  contour  confidence  regions  were  generated  for  each 
target's  least  squares  solution.   Using  a  confidence  level 
of  90  percent,  national  network  confidence  regions  contained 
90.38  percent  of  the  generated  targets.   Results  for  the 
global  network  were  90.44  percent.   Such  close  agreement 
between  the  confidence  level  and  containment  in  the  simula- 
tions is  dependent  upon  a  reliable  standard  normal  azimuth 
error  generator.   In  production  use,  such  reliability  could 
only  be  achieved  by  having  a  good  estimate  of  DF  site 
variances. 

A.   NATIONAL  NETWORK 

Figure  26  is  a  scatter  plot  of  confidence  region  area 
versus  location  error.   Again  using  a  common  scale  to  aid 
network  comparison,  eleven  points  off  the  scale  are  entered 
above  or  to  the  right  of  an  indicator.   Although  there  is 
one  hopeful  cluster  near  the  plot  origin,  the  large  variance 
in  location  errors  at  given  confidence  region  areas  is 
indicative  of  the  futility  of  predicting  location  error 
based  on  the  size  of  this  reliable  confidence  region  model. 

Figure  27  is  a  multi-level  histogram  of  location  error 
at  the  various  degrees  of  freedom,  demonstrating  the  same 
sort  of  correlation  evident  in  Figure  13.   Figure  28 
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illustrates  the  correlation  between  location  error  and  the 
mean  acute  angle  of  intersection,  present  even  with  mixed 
degrees  of  freedom.   The  same  holds  true  in  Figure  29,  a 
plot  of  location  error  versus  the  mean  arc  statistic. 

B.   GLOBAL  NETWORK 

Figures  30  through  33  are  similar  plots  for  the  global 
network.   On  Figure  30,  location  error  and  confidence  region 
area  appear  independent.   Figure  31  shows  good  correlation 
between  location  error  and  degrees  of  freedom.   Figures  32 
and  33  show  nearly  complete  masking  of  the  previously 
established  correlations  due  to  the  layout  of  the  network. 


80 


■      I-     -1      r- 

' 

"T 

X   -L    ^ 

j_ 

1 

i 

1 

i 

!      1 

i 

i 

i 

! 

1 

IM 

i 

t: 

j 

1 

1           Inn 

1 

-j*-^ 

I  ! 

1 

i . 

-                               ! 

• 

| 

X 

T    ■ 

1 

1 

I 

1 

L_ 

1 

1 

i 

mX        i 

ii 

. 

.  i 

1  f    i 

,  . 

i 

I  I      i  -><-n 

i 

i    . 

iCO-'C^CP^    p.<..o  f 

,\  erf.   Lo<uuL.\.*jy 

n     QJVlQ'L 

i 

JJtU 

I    , 

a 

-il^LC  ^4*"  * 

■Lo-i. 

•i-C,    AJl 

_L_. 

1 

lat-lor 

xz2L/ie 

'Mia-^L  . 

1        1 

1    1 

1 

1 

M-      -4-X 

1 

It 

| 

| 

1     [ 

+;  • 

1                                                   1 

JuO 

i 

1 

i 

1 

n 

i 

| 

i 

1        . 

1 

|                 1 

|__ 

II 

~r 

I       xx 

i 

r  i 

I 

r 

1 

or/i 

|  J.JU 

j 

1     1     1 

| 

' 

M 

i 

!       ! 

i 

1    1 

X 

1  1 

_ ,. 

I 

4_ 

1 

i  i 

1 

1 

i    1 

1            '  on/i 

t          '  -r^" 

mi 

i 

i        i               Mm 

. 

i 

II     i  ' 

1  i 

1 

l 

f  I  T      T 

1       II 

i   i 

j_ 

1  i 

|       :       1       , 

ill!     i  i 

1 

• 

i 

1 

i     !     II     MM  1        i     1 

i 

X 

i  i  • 

!    " 

1                  1  i  1          i     1 

i               1 

1 

i  | 

* 

1    \*'    • 

|     1     1        |  1  |  1             1 

1                 !    ! 

.        I|        1 

1    I       1   1                 Mi 

150 

i  ;  ;  ' 

1  '  1  1 

i      i      .  '->u 

! 

i 

i 

•    ; 

±±          T 

..  it. 

1 

II 

H         * 

iii     i  i  i  i     1              i  i  i 

ut    r 

4- 

1     1 

i  i 

1  i 

1 

l 

r 

Ml       Ml          ;   i                Ml 

x 

i  | 

U 

mi                 i  |          |     M 

1 

m 

1    1    1 

• 

i         .    i 

i 

I'M 

j  |  i 

I 

• 

!    i 

1       II            1  i  1     1  1 

i  • 

i  | 

hi 

IT 

-• 
1 

i  •• 

ill       1                  MM       1       1 

i 

1   :p 

i 

I 

1 

I  1 

ill 

ill       I                     MlM 

I  I  I  I     ' inn 

ii' 

:   ! 

1   !   i   i 

I'll     ! 

\             tuu 

• 

i    LI    j 

• 

1             1     i 

i   ! 

X 

1 

[ 

• 

ii  j 

tl     1                       M     j     i 

! 

; 

1 

T> 

i 

11      1  ' 

•  • 

IM           |                         1       1       I 

M-          ±^t 

Mil 

•            i   «, 

Ml 

• 

||           1                  II1 

1             ! 

'ill 

1  1  i 

.. . 

. 

i  *  1 

T 

_m 

j 

1 

1 

i 

T  i  ! 

L_ 

i  1 

!     i 

1   1 

7 

II            1 

11             ! 

MM 

1 

1 

i  l 

.vx: 

4~- 

•  i 

-| 

1" 

ill            ' 

1  1  i  1 

i 

1 

! 

1 

i 

-+*■    -- 

• 

MM              i  i 

'  rn 

1  '  '<  • 

• 

". 

\  . 

i_ ,__     50 

1 

j    i 

i«i  ,  • 

•»  ♦. 

■ 

1  i     i 

Mil 

1 

1    ! 

1    1 

• 

.  -,i 

*l 

* 

•     |  i  " 

II             1  M     i 

i 

i   i   1 

1 

i 

1    i 

•  i 

i  i  :  : 

1  i  i  • 

'                  1    !        : 

! 

, 

i 

L 

'  1  '•'■ 

i  i  . 

:* 

••• 

•      • 

III        ||                   Ml        1    1    i    1 

! 

1  i  i  i 

•        * 

••    •. 

*•• 

i  !  i 

1  f  1 

•  • 

* 

• 
•  •• 

•*• 

i  1 

!   I   i 

i  '  ! 

*•  . 

.         • 

♦i 

.  • 

1  III 

1      :      M      ;      I'M      ' 

'                • 

! 

• 

!*• 

•  - 

.     '  ,,  ' 

i    j  1    1    1    1  1 

1  ! 

•    *H  J 

•    . 

•  ■ 

• 

Y) 

•  * 

• 

•  •  ■ 

!  ;             ■ 

2 

i 

At 

> 

.  &it 

i 

i-Jlii 

1       i 

4 

' 

HK 

' 

i  !     bT 

Clh: 

i/cjt''e&a 

!    i    ii 

■ 

III 

■    i 

II     ii.  ''  '      Mil      1  M 

II 

1 

i  . 1    

Z1QURL 

29  \ 

1    1    ! 

-  -~ 

..... 

Mil 

|_    .j.       1           |        |  I      ■ 

ii         i 

1 

1    <    ^_ 

1  i 

1  1  1  i 

i  1 

1             1 

i    X 

_i          i        M 

'    ' 

1 

II 

i 

1    1 

i        i  i    i  i  i  i    i 

1           ■ : 

1  1 

I  1  ^ 

1    1    1    1 

1   1       ! 

;  ! 

-          T  _r              ~       ~     -        T              I         -  -F-F 

1 

Zp-IL"1        -J j _                             _L                                          !  Tj                         _L                                   J_     It"       ~     "'     ~~ 

!           T  i  -  !  t ..  .                       _  .         .  _b     : 

/ff/c L 

1     '     '^nn                                           '     ' 

<tUU                                                                                                                                                                                                                 —  --■^- 

1 

i                                                                                                                                     1 

L      _L           i 

1  i LL                                                  .                              J   .     L. 

iirn                                                                               Loca^Xott  eAX,o4.1  </&L<>i.4 

—  acyy-tee^    o^,   {si<z.&cLoifi-  ^rx 

ri IjnJbcL^    no <t±i  ci/i k~ 

.... ._.    .                                                                          ......    _L       . 

I                  :  *>nn           I                                                                               III! 

\                                                                                                                                                    I 

| 

i    iy               _                                    1      '      " 

I 

i 

i 

i  i    _j_ 

\           6 nn        I                                                                                                              j 

•ou 

i 

I                      ?,w.        :                                                                     II                                            III!                            I    ,                !                        I 

i            .iuu 

ill                                                   |     i                                                          I        i     i  i             I 

I                                                                                                                          I        i     |  i          I  i 

i                                          !                   L 

i      i 

l                          III                                                                         Li                                               :     !     I                     |         j                      I    |                 '     j    I             I    l    '    I 

i                                                    ii                           I  I               !     i                         III 

II                                                 M                                  I                         I                    i  i          i  ■■  i 

I                                I                 I  i            i                    i  I            ;  i            i                    i            I  I 

1rri                  '                                                                                                                                                                                                    I         I    |         I         I     I 

I    I    i               ' JW                                                                                                                                                                                                                     •        i        i    ' 

II                        I                                                                                 I        i        I    I            i                ||                    I        I            I            I    |                        I        L 

i                                         I                                '        I            I 

I      i      I                                            ™,                                     I   i                      !   L   i             I          !                I      I             i      I          |   M 

~l                                                                                                                                     I         I      I 

I    .  j                                                                                                      'I'll 

mr\ 

0         I                                                                                                                I      I      i  l  i           l  i 

M                                              i        I     I     L                            i                       I                         I  I 

1   I     I             I      •            I        l        i                    I                  I        i  I  i                               I               I 

l                    I        i                                                             III 

>        I           |            [    |    I        II                   II               I                                                         l        i        I    j    I            !        I 

i        I    I                           i                         .III                                                 l|L            |        I 

I  I     I     I       I       mi       M              !                              III 

Hi       I     I     I                II     i  I  |  I     I  I  i  i         III    mI     I       i 

: rn                                                                  II                                                                 I      i  I  I  I      ill!      '   :   '   i 

J                                                                                   I                     I                             L  *  .  t   I    i     ._._  J                             |    I    i    l       I    i 

I                                                                 |ll                                                         l|l||!. 

j                                                              ,_[                     I    I    I _..  _U                                                                                              ■!       I    |       I    !    I 

|  i               i                                I                                                                 I      i  »      till 

I  111        1    1                                            n 

II                                    1                                    1    !    |        i    i    1            |    1 

1                                                                                                                                ' 

1    1                    1 

L/i                                                ~l 

\Y                              1           \i   \         .  <;_.          f, .          "'             8           o  .                            ^i/ 

!  I  !        I                                    •'.  i  !  I    *'  I  i  !  i-'  i  i  i  I  *>  i     |  \t        i     '        I  ii         |                    A^f     _U 

I   |                                                    i          )                       I   I      I   I  I                I                                 I                I             I             |.  i 

:                                                                                          I        I        r^/iir      o,                   '    ]                                                  i- 

TO$ll£L  11 

II                                                                                                I                       I              I    j         I         I     I         I         |             j             |         I                                | 

II                                     I  I                     i       i  I     I       I     I            I'll                       I  I j  -     . 

i     — ' — I — I —  t  I    I  i                   t    J    h                I      I  ~r~\~          ._.    it  ! 

,         I     I     I     !  I     I  i  !           i           i              I     i  i     i                                                 LL  i       i 

83 


ii±#nr~=~: 

! 

I 

"    i 

L 

— + 

:=£:££=- 

.m        _l  . 

_L    . 

T 

ii 

1 

1      +          r 

.  Jl\fSL        It 

i 

|       I 

41 r          | 

1 

| 

1 

4on          ]        ' 

i    j       i 

i  HUU 

i      ill            |  i 

l  i      I 

1       1 

1              11 

I 

T^    | 

It  It 

t      _t 

:  ztzrr 

i                 i 

r         i            i  i 

IE    f 

.A    _1_        It 

I 

i 

.     .      ^L 

it           niz; 

, I.,  i 

i 

■yrn        Ml                   1    '    be 

aZtaz   px.o  t  of.  r^.occw^can. 

ejiA,on, 

i  J  JV 

Jt>su<-m.e.an  acutQ-a/i 

nf    .         |    1 

Of- 

i   in 

EkJu  e.c.Lin  i  An   a  LahaL  n 

e  iuo  -tk. 

II           ill      i 

!   i      I   i 



1                i 

li      I 

1 

1      1  1 

.' 

1       _L 

4     .  ■ 

1 

1  inn     '    '  ' 

1  i 

i      1    " 

| 

JUU 

! 

1 

I  j 

1    ! 

.J    1 

|    1 

1    1 

1 

1 

.    _..  it 

i 

i    1 

1 

: 

1 

it         ' 

~4— 

j    4- 

1    '    i 

'J'Tfi 

ffffl 

! 

1 

1                | 

i 

4-     T- 

1    '■    \ 

! 

-     -     it     4- 

1  i  1 

i 

■  +       ^4^ 

i 

1 

^~ 

i  i 

i 

II 

1 

I  I  , 

1 

rn 

]      1 

T  - 

J    j_ 

TI 

1 

i 

. 

i  i 

Inn         l           ! 

i  1          1  |  i 

1 

r  r  1  4- 

•CUU 

i     i 

_  1  I 

!        1 

:    .    ±    t 

r 

zu 

1 

i 

i 

in 

~r 

1 

1      j       1 

i        i        l 

i 

1 

:      '  in 

1 

i    |  I 

|    1    |    !                               !                i 

i    M    i 

[ 

i 



1  i 

1    1 

!  i                    1 

i 

i  1 

1 

:     I 

1                      i 

i 

1 1 

i 

|l                                   -J uLi 

I 

1 1 

fpr,         1     i         !         j 

i  | 

I^U                  ill             li 

i 

1  i  1  1 

i        1                III 

i 

i 

11 

| 

i 

1    !     ! 

i 

i 

ii 

'  I 

-r4- 

1     1 

! 

i 

i         I  i 

i    i 

1  1 

1 

[_ 

1    '    !    >                                            1 

1 

i  i 

1 

!  ; 

1 

i               ■                       '                                      ;                       ! 

1    •                         1          1     1 

!      I 

Ml; 

;   l                   1 

1               III 

• 

it         ill 

i 

i   1   1 

i 

i                       ; 

J 

i 

__u 

1 

III                       j 

| 

I         L 

1 

1 

1  1 

inn 

1 

1  !            ! 

1 

100            :      |      1      :            ill            1 

!               |     ' 

II 

ill 

1  1 

■   i   '      l   1   I   i      1         !      ! 

!    1 

I 

1              II 

!  !     i     1             -i 

• 

1 

1    1 

T 

\         1         1 

III 

1             •     1 

. 

1    1           i 

1  1 

II      1    i 

1               1               II 

1 

I         1         1 

1    1          1        1  1 

1       1           „l.|    I'1 

•  -        ,  . 

I         :         1    ! 

|         | 

1                 1           Ill 

"j  i    mM 

'|*       :     |           1 

1    | 

T 

1 

L 

.     !    1     •!  "    1 

. 

!  1 

r     , 

1 

: 

IE 

II        1            1 

Ltl 

I-1     4       •  I 

1  i       1 

IT 

\    \rn             \    1             1    1 

... 

I      ,    |,    «, 

li     j  •  i 

i 

1   .50 

♦                                             •        •*•     *^a             i. 

;    : 

1 

1 

i     1                <           ♦                       i 

'     i     |           :     ;     ! 

1   !  I 

11     11 

i 

*• 

* 

i       ! 

±-J^ 

1 

1       !    1   _ 

1 

III 

'   . 

;J__L 

1           1 

1 

j_ 

Mill           1 

*  •.  ;■ 

*•♦« 

i           i 

!   !      1 

_!3' 

• 

'-?*.• 

1         i  I 

•,  rrv 

r^ 

r*     , 

>  < 

• 

i     i 

! 

, 

1   i                                        JL 

.  *  ' 

.■£,« 

%' 

> 

i 

|     i           i     i     | 

1 

1  I 

i" 

•:. 

i 

1  1 

1  . .  1 1 

i 

'  n 

•  t     ♦» 

\  \  P                          ">? 

•  i      !      4i>- 

fin 

f 

>££> 

:ce^3 

I      i   i 

u.y 

i  r; 

i 

i 

t 

1 

!      !   1 

1  '     i 

w 

;m 

R 

1 

i      1   '   ,   . 

1      !  ! 

73£ 

2 

~        1 

| 

I 

1   I  1 

i       !   1 

Mill 

I               i  1          i| 

I  1  1  I 

4 

:    4  ,4! 

1 

1 L4         - 

i 

i  i     1 

_           i III 

1  . . 

I 

4~  '  - 

i               i   i   !         | 

I               :   ! 

1   II 

.           1 

1    Mi- 

1      1  i 

-14- 

_L 

....       , 

i   I            i 

84 


I_                             ■"" 

■ 

i_         in 

! 

-     -              _|_        X        J 

'  X      x"          bh 

x     ~  ±            iX 

H — LZ     __                    . t  x       : 

X 

nm iii-    -                           .1 

run 

x 

•inn     I                                     ^     ' 

i 

x  : 

1 

t    : 

-                                            M           Tl 

1 

-              ±„riiii^- 

j  *  j 

. 

i     I 

ft                  _____     ______  "    r  ' 

It^                            1  .iou.tte^i  p.Lol 

3/:   -cocax^; 

m.   ttsipscv 

JjU                                    .  ,      M 

IG    -O-i 

ta^co- 

Ud  ■ 

* 

L/T 

...  .  -                                    tJaJSKL  ^nc-two-t 

t-     - 

r               i 

l 

1 

I 

'  inn                                                 ' 

JUU                                             i 

X               -r11- 

-   .                 ..    _,                       .             _i         __   . 

!            i 

i 

1 

OM             ■                      ' 

I  i     i 

Iju                     I 

1       1            i                '     '       1 

1 

!                                     1  \ 

i 

1                    |      1  l  |        ■      1 

1 

III                 1 

1             1                      I      M             |    H^ 

1                                     II             i      1 

nnn 

U.UU 

i  |     Ml            1  i 

II        i  1          i        i 

Mill           i 

i      :m             i  i 

1 

1      II      i        i      1  M        1  I  i 

1          1          I              MM 

II          Ili                                         I 

1 

|                       '     |       !     '  : 

jrrf       \    \                  M               i 

i 

7-^               I                                  !               | 

i 

1 

i 

i            1  1               i  i   '      II   II      M 

1 

1                             ''ill 

1                                  1          I                                                 !     ! 

1                                                                            -                               II 

i  i 

l            i    j    :        i    ;    i    .        Ill' 

i 

1              1     j                  !     1     1 

II         III              i     |     1     i 

| 

j                                                       !|         Ili              1     i     1     i 

mn 

IUU          | 

I                          I 

M      till        |  1 

■ 

.  *  i 

1       M       II       i   1   |   1 

« 

. 

ill 

i 

H       1 

Mil               1       II 

||                  I                      ! 

t 

4 

>  i  • 

• 

-"  i  j— : 

»        . 

I               i   ;                     j      1                 .XI 

-M-H-' 

► 

l               1 

I 

" 

•• 

rn 

. 

•*       J 

ju            i                 i              I,  |     | 

»j    ••• 

..    t. 

.  •     •                     I   M         |  i 

i          i                        I  '     i 

__ 
* 

M  *  i* 

*•    i  * 

i                          !     i 

»  m 

^       l« 

•    • 

i     1     1             III                 1             III'         1 

,\ 

•    -  •  —  • 

♦  _.  •*• 

»M                 II            1              i   '      '   ; 

I 

•    •! 

■ 

.'  i  !  V'< 

^:.#; 

•     ♦•           #^ 

♦' 

■  «      * 

-  •    •• 

ur 

;»i 

I                           I              !              ||              I                                     III 

>  1  ••*! 

ii          ill     II  '                   .     - 

!         i         .MM                   i  II  I      I 

_'! 

•  •  • 

►  .    •     • 

II          III       i 

Ia                                             

• 

Dl                        LCl 

N 

1 

J.)                        **\) 

i    |    1 

l   °< 

1          I    Ol/                           1     i/K-J.'  _{_,t_-0 

T1 

mk 

->-> 

•  1  ! 

.^_          I          

f-J 

33 

ill! 

I        1            I    1                III! 

MM         I  i  |      j 

!                         1     1     1        '                 1            ___    |.J     . 

1         I      |                        1      !                  !   1   ! 

i 

I       I         L          in  '  -      - 

i             x  m              !  .  - 

X         i  M i.l_ 

~T 

i  i 

1     | 

85 


IX.   CONCLUSION 

The  Fix  Program  and  Chi-Square  contour  confidence  region 
subroutine  could  be  integrated  into  the  software  of  position 
location  systems  by  direction  finding.   The  present  program 
(with  target  and  azimuth  error  generators  included)  could  be 
used  to  simulate  an  existing  DF  network  for  accuracy  studies 
on  potential  targets  in  different  regions  of  the  globe  or  for 
studies  on  optimal  placement  of  DF  sites  in  a  proposed  net- 
work. 

The  relationships  demonstrated  in  the  sections  on  loca- 
tion error  analysis  indicate  that  larger  location  errors 
accompany  contraction  of  the  network  layout;  however,  this 
factor  would  be  mitigated  by  a  practice  of  assigning  regional 
rather  than  globally  distributed  targets  to  the  network  for 
prosecution. 

Countering  the  widespread  belief  that  small  confidence 
regions  indicate  accurate  locations  and  large  confidence 
regions  indicate  poor  locations,  it  has  been  demonstrated 
that  for  a  reliable  model  location  error  and  confidence 
region  area  are  nearly  independent. 

The  relationships  examined  as  parts  of  location  error 
analysis  suggest  an  alternative  answer  for  the  question  of 
how  accurate  is  a  given  location.   For  an  existing  network, 
data  could  be  collected  from  locations  of  known  transmitters. 
Such  data  would  include  statistics  on  number  of  signal 
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azimuths  available  to  make  the  fix,  mean  acute  angle  of  inter- 
section and  mean  arc  statistic.   With  known  location  error 
as  the  dependent  variable,  multiple  linear  regression 
analysis  could  be  used  to  develop  regression  and  correlation 
coefficients  for  predicting  location  error  as  a  function  of 
independent  variables  measured  or  calculated  at  the  time  of 
location . 

The  workability  of  such  an  approach  was  not  demonstrated 
in  the  present  paper  for  several  reasons:   (1)  The  effort 
required  to  perform  the  multiple  linear  regression  (using  a 
standard  program  from  the  IBM  scientific  subroutine  package) 
would  be  a  project  in  itself  and  would  be  best  reserved  for 
actual  data  from  an  existing  network.   (2)  System  operators 
would  be  aware  of  other  measurable  factors  for  inclusion  as 
independent  variables  in  the  regression  analysis,  such  as 
variance  on  given  signal  azimuth  measurements,  signal 
strength  and  measurement  duration. 
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FIX  PROGRAM 


0900 


9 
1 


RFAL  LAMBD 
D I  MENS  I  ON 
STAA7( 19) 
COMMON/ SIT 
COMMON/ sir, 
DATA  TWOPI 
DEGR AD/1 .7 
PIHALF/1.5 
CATA  DFSIT 


IDENT 
1., 
2., 

3.  , 

4., 

5., 

6., 

7.* 

8., 

9., 

10., 

11.  , 

12., 

IB., 

14.  , 

15., 

16.  , 

17.  , 
18., 
19., 

READ 
READ 


LONG 
210. 
173. 
288. 
262. 
254. 
276. 
253. 
247. 
202. 
241. 
279. 
182. 
266. 
273. 
269. 
283. 
240. 
2  37. 
262. 


1000 
C  PERF 
C  AND 


{ 
(5,1) 

IF  (NUMBER 
WRITE  (6,9 
FORMAT  CI 
FORMAT  (12 
DO  1000  K 
READ  (5,1 ) 
THETA(K) =T 
S  I  T  E  (  1 ,  K  )  = 
SITE(2  «K)  = 
SITE( 3,K)  = 
SITE(4,K)= 
SITE(5.K) = 
SITE( 6,K)  = 
S I  T  E  (  7  ,  K  )  = 
SITF(8,K)= 
SITE(9,K)= 
SITE(10,K) 


A,N, 
DFSI 

ES/S 
NAL/ 
/6.2 
4532 
7079 
E/ 

117, 
250, 
933, 
600, 
983, 
950, 
517, 
983, 
133, 
750. 
800, 
633, 
783, 
200, 
383, 
767, 
200, 
667, 
667, 
KIL 
NUM 
.EG. 
) 

'  ) 

,F8. 
=  1, 
ID, 
HETA 
DFSI 
SINK 

cost 

AT  AN 
SINK 
COS( 
DFSI 
DFSI 
DFSI 
=  DFS 
ATIiM 


NL,NP,NLL,NPP,NPL 

TE(6,19)  ,D0MEGA(19) , SYMBOL ( 4 ) , BET  A (  19), 


ITE( 
THET 
8318 
9E-2 
6/,SYfBCL/«N' 


NUMBER 


10,19) 

A  ( 1  9  ) 

5/,PI/3. 141593/,RADDEG/57.29  578/, 

/  ,  FNC ECC/ 1. 00 676 S/, SCREEN/. 9999999/, 


W  / 


LA 
61 
52 
42 
27 
39 
42 
31 
46 
21 
34 
25 
28 
44 
36 
29 
33 
39 
47 
37 
CU 


T 

.217, 

.533, 

.350, 

.  ECO, 

.717, 

.333, 

.750, 

.600, 

.217, 

.050, 

.767, 

.  217, 

.983, 

.150, 

.967, 

.667, 

.  517  , 

.600, 

.633, 

T 


SITE  NAME 


'ANCH 

, 

OR  AG' 

,%  E    '  , 

'ATTU 

, 

' BOST' 

, 

ON   « 

•CORP 

T 

'US  c 

1'HRI S«  \ 

•  DENV 

» 

»ER   •  , 

•DETR" 

, 

CIT  ' 

'EL  P 

, 

'A  SO  « 

•HELE 

, 

1 NA   ', 

«HCNC< 

, 

LULU* 

'LOS 

, 

•ANGE' 

!  J  L  E  S  • ! 

' MI  AM" 

, 

I    •  , 

•MIDW 

, 

AY   ' 

•MINN 

, 

'EAPO' 

! 'lis  • ! 

•  NASH' 

, 

VILL' 

r  »  p    1  % 

•  N  E  U 

, 

!ORLE» 

t «ANS  ' , 

•NORF 

, 

'OLK  ' 

•RENO' 

, 

'SEAT 

, 

'TLE  ' 

•WICH' 

, 

ITA  ' 

!'   •/ 

BER 

0)  GO  TO  4200 


3) 

NUMB 
THET 
(  K  )  * 
TE(2 
SITE 
SITE 
(0.9 
SITE 
SITE 
TE(1 
TE(4 
TE  (  5 
ITE( 
G  GO 


FR 

A(K) 

DEGRAD 

,  ID)*DEGRAD 

(1,K)  ) 

(  1  ,  K  )  ) 

5  327  73*TAN(DFSITE(3,ID)*DEGRAD) ) 

(4,K)  ) 

(  '. ,  K  )  ) 

,  ID) 

,  in) 

,  ID) 

6, ID) 

LDEN  SECTION  SEARCHES  ON  LONGITUDE 


1100 


1200 


ORM  ALTERN 

LATITUDE 

PSI=0. 

LAMBDA =0 

SAVEL=l AM 3 

SAVEP=PSI 

DEGF=61560 

CO  1400  J 

DIV  =  J 

X=1./DIV*TW0PI 

Y  =  0. 

CALL  SEARC 

IF  (KILL.E 

CALL  SHIFT 

CONTINUE 

x  =  o. 

Y=1./DIV*P 

CALL  SEARC 

IF  (KILL 

CALL  SHIFT 

GO  TO  1300 

CALL  AZDISTt  SA VE L . SA VE P , LAMBDA , P SI  ,AZGO, CIST) 


IDA 

0. 
=  1, 


H(LA 

1.2) 

(  SAV 


I 

H(LA 
Q.2) 
(SAV 


MBCA 

GO 
EL,S 


MBCA 

GO 
EL,S 


,PSI,X,Y,DEGF) 

TO  1100 

AVEP,LAMBDA,PSI ,1,DEGF,DIST) 


,PS I,X,Y,DEGF) 

TC  1200 

AVEP,LAMBDA,PSI, 1,DEGF,DIST) 


88 


1300 
1400 
1500 
160  0 


C  EVA 


1700 
1803 
1900 

2000 


2100 
2200 

2300 
2400 


SAVEL=L 
SAVFP=P 
IF  (OIS 

CCNTINU 
M  =  0 
M=M+1 
COSLAM= 

SINLAM= 

COSPSI^ 

SINPSI= 
FL  =  0. 
FP^O. 
FLL-=0. 
FPP=0. 
FPL=0. 
LUATE    FI 
DO    2450 
A=l. 
C=0. 
UA=SINP 
UB=COSP 
UC=CQSL 

un=-siN 

UE=COSP 

UF=SINP 

UG=UF+U 

UH=UB-U 

UI=-UE* 

N=UA-UP 

NL=-U3- 

NP  =  UE-«-U 

NLL=UB* 

NPP=-N 

NPL=UF* 

IF     (UG. 

IF     (UG. 

0=SIN( A 

UJ=1./D 

DL=UJ*U 

DP=-UJ* 

DLL=UJ* 

DPP=UJ* 

DPL=UJ* 

COSB=N* 

IF    (COS 

IF     (COS 

E=ARCOS 

COTB=CO 

CSCB=1. 

BL=UJ*( 

BP=UJ*( 

UK=UJ*C 

BLL=UK* 

1BL*NL*C 
BPP=UK* 

1BP*NP*C 
BPL=UK* 

1 BP*NL*C 
TEST=SI 
IF  (TES 
IF  (TES 
IF  (TES 
A=-l. 
C-TWGPI 
TEST=A* 
BETA(K  ) 
IF  (ABS 
IF  (TES 
C=C+TWO 
GO  TO  2 
C=C-TWO 
CNEGA=A 


AMBDA 
SI 

T.LE.500. 
E 


CCSCLAMBC 

SIN(LAMBD 

cnstpsi ) 

SINCPSI  ) 


RST    AND    S 
K    =    1,NU 


SI*SITE(6 

SI*SITE  (5 

AM* SITE  (3 

LAM*SITE( 

SI*SITE(6 

SI*SITE  (5 

E*UC 

A*UC 

UD 

*uc 

UD 

F*UC 

UC 

UD 

GE.l.)    UG 

LE.-l.J     U 

RCOS(UG)  } 


)    GO    TO    1500 


A) 
A) 


ECOND  PARTIALS  OF  OBJECTIVE  FUNCTION 
MBER 


,K) 
K) 


,K)+SINLAM*SITE(2*K) 

3,K)+C0SLAM*SITE(2»K) 

,K) 

,K) 


=SCREEN 
G=-SCREEN 


G#UI 

UG*UH 

(UG*U 

(UG*U 

(UG*U 

UJ 

B.GE. 

B.LE. 

(COSB 

TAN(B 

/  S  I N  ( 

DL^CO 

DP*CO 

SCB 

(UJ*D 

OTB) 

(UJ*D 

OTB) 

(UJ*D 

OTB) 

TF(  1  , 

T)  17 

T  +  PI  ) 

T-PI  ) 


E*UC 

G-UH 
A*UD 

1.  ) 
-1.  ) 

) 
) 

B) 

TB-N 

TB-N 

L*  (  N 

P*  ( N 

P#(N 

K)-L 

00,2 

190 

190 


+UI*UI-DL*DL) 

*UH-DP*DP) 

+UH*UI-DP*DL) 

COSB=SCREEN 
COSB=-SCREEN 


L^CSCB) 
P*CSCB) 

L-DL*COSB)+DLL*COSB-BL*DL*CSCS-NLL+ 

P-DP*COSB)+CPP*COSB-BP*DP*CSCB-NPP+ 

L-DL*COSB)+DPL*COSB-BP*DL*CSCB-NPL+ 

AMBDA 
000,  1800 
0,2000,2000 
C,2CCO,2000 


B+C-THETA 

=  (A*B+C  )* 

(TEST)-PI 

T)  2200,2 

PI 

400 

PI 

*e  +  C-THETA(K  ) 


(K) 

RACCEG 

)  2400,2400,2100 

400,2300 
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2450 


C  IF  E 
C  NM  n 
C  VECT 


2500 


2600 
C  REFI 


2700 


C  CHEC 
2800 
C  S  FAR 
C  IS  R 


2900 
C  IF  T 
C  SECT 


3000 
310C 


3200 
3300 

3400 
3500 


DCMEGA(K)=OMEG 
UL=2.*A#(A*B+C 

FL=FL+UL*BL 

FP  =  FP+UL=?B0 

FLL=FLL+2.*BL* 

FPP=FPP+2."RP* 

FPL=FPL+2.^BP* 

CONTINUE 

UM=l./( FLL*FPP 

X=(FPL*FP-FPP* 

Y=(FPL*FL-FLL* 

TTHER  X  OR  Y  A 

F  THE  INDICATE 

OR 

IF  (AMAXK  X,Y) 

G0X=LAM3CA+X 

GCY=PSI+Y 

CALL  AZDIST(LA 

IF  (DIST.LE.30 

CONTINUE 

DIST=SORT(X*X+ 

X=X/DIST 

Y=Y/DI ST 

CONT  INUE 

NE  INDICATED  M 

CALL  SEARCH(LA 

IF  (KILL.E0.2) 

CALL  SHIFKSAV 

GO  TO  2  9  00 

CALL    AZDISKSA 

SAVEL=LAMBDA 

SAVEP=PSI 

K    FOR    MOVE    LES 

IF     (DIST.LE.CU 

CHING    STOPS    AF 

ECU  EST  ED 

IF     (KILL. EG. 2) 

IF     (M.E0.30)    G 

CONTINUE 

HE    SOLUTION    HA 

ICN    SEARCH    AFT 

IF     (M0D(F,10). 

X=PIHALF 

Y  =  0. 

CALL    SFAPCH(LA 

IF    (KILL. EG. 2) 

CALL  SHIFTtSAV 

GO  TO  3100 

SAVEL=LAMBDA 

SAVEP=PSI 

X  =  0. 

Y=PIHALF 

CALL  SEARCHtLA 

IF  (KILL. EG. 2) 

CALL  SHIFT(SAV 

GO  TO  1600 

SAVEL=LAMBQA 

SAVEP=PSI 

GO  TO  1600 

DEGLAM=LAMBDA« 

IF  (DEGLAM.GT. 

X=SYMB0L(3  ) 

GO  TO  3500 

DEGLAM=360.-DE 

X=SYMB0L(4) 

NDEGX=DEGLAM 

FDEG=NDEGX 

NMINX=( DEGLAM- 

DEGPHI=ATAN(FN 

IF  (DEGPHI.LT. 

Y=SYMBOL( 1) 

GO  TO  3700 


A*RABDEG 
-THETA(K)  ) 


3L+UL*BLL 

BP+UL*EPP 
BL+UL*8PL 

-FPL*FPL) 

FL )*UM 

CP  )*UM 

RE  GREATER  THAN  PI  AND  UNTIL  WITHIN  30 

0  MOVE,  CONVERT  (X,YJ  TO  A  RADIAN  UNIT 

.GE.PI)  GO  TO  2500 


MB DA , PS  I  ,  GOX,GOY,  AZGCDIST) 
.  )  GO  TO  2600 

Y*Y) 


OVE  BY  A  GOLDEN  SECTION  SEARCH 
MBDA,PSI  ,X,Y,DEGF) 

GO  TO  2700 
EL, S A VEP, LAMBDA, PS  I, 2, DEGF, D I  ST) 

VEL, SAVEP,LAMBDA,PSI,AZGO,DI ST) 


S  THAN  MINIMUM  STEP 

T)  GO  TO  3300 

TER  30  ITERATIONS  UNLESS  SEARCH  PROGRESS 

GC  TO  2900 
0  TO  3300 

S  NOT  BEEN  FOUND,  APPLY  THE  GCLDEN 
ER  EVERY  TENTH  ITERATION. 
NE.OJ  GO  TO  1600 


MBCA.PST ,X,Y,DEGF) 

GO  TO  3000 
EL,  S A VEP, LAMBDA, PS  I, 1,DEGF,DIST) 


MBCA,PSI,X,Y,DEGF) 

GC    TC    3200 
FL,S A VEP, LAMBDA, PS  I , 1,DEGF,DIST) 


RADDEG 

180. )    GO    TO    3400 


GLAM 


FCEG)*60.+.5 
CECC*TAN(PSI ) )*RADDEG 
0.  )     GC    TO    3600 
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) 


4CJX 


•  S 

A 
30 
NO 


360C       Y=SYMB0L(2) 

DEGPHI  =■— DEGPHI 

3700       NDEGY=DEGPHI 
FDEG=NDEGY 
NMINY=( DEGPHI-FDEG 
WRITE     (6,2)     NPFGX, 

2  FORMAT     ( «0' ,30X,I3 
1'     LFAST    SQUARES     SO 

WRITE     (6,3)     DEGF 

3  FORMAT     (31X 
1'     SQUARED.' 

WRITF     (6,4) 

4  FORMAT     (33X 
WRITE     (6,5) 

5  FORMAT     CO' 
1  «     ID    INTERCEPT 
2*     THETA    VALS',/ 
3'=    OMEGA       NM 

SUMARC=0. 

DC    4000    K    =    If  NUMB 
ID=SITE( 7,K)+. 5 
DTFETA=THETA(K )*RA 
CALL    AZDISTtLAMBBA 

IDT  ST) 
NMSITE=DIST+.5 
ARCC=(DIST/3439.92 
ANGLEB=ABS(DOMEGA( 
ARCB=ARSIN(SIN( ANG 
ANGLEA=ARCOS(TAN(A 
STAT=ARCC*RADDEG 
IF     (STAT.GT.90. )     S 
IF     (DOMEGA(K) .GT.O 
AZTHET=SITEAZ-ANGL 
IF     ( AZTHET.LT.O. ) 
STAAZ(K) =AZTHET+9D 
IF     (STAAZ(K).GT.36 
GC    TO    3900 
AZTHET=SITEAZ+ANGL 
IF     (A7ThET.GT.360. 
STAAZ(K)=AZTHET-90 
IF     ( STAAZ(K) .LT.O. 
NMTHET=ARCB*34  39.9 
WRITE     (6,6)     ID,  SIT 

1DTHETA,DCMFGA(K) ,N 
FORMAT  (30X.I3.1X, 
SUMARC=SUMARC+ STAT 
DNUM=NUMBER 
AVGARC=SUMARC/DNUM 
WRITE     (6,7)     AVGARC 

7  FORMAT     (  «0« ,30X, »M 
CCUNT=0. 
SUMANG=0. 
NEND=NUMBER-1 

DO  4100  1=1, NEND 
JBEGIN=I+1 
DO  4100  J  =  JBEGIN 
ACUTE=ABS(STAAZ(I ) 
IF  (ACUTE. GT. 180. ) 
IF  (ACUTE. GT. 90. ) 
SUMANG=SUMANG+ACUT 
4100       CCUNT=CGUNT+1. 

AVGANG=SUMANG/CCUN 
WRITF     (6,8)     AVGANG 

8  FORMAT     (31X,'MEAN 
1»     DEGREES.'  ) 

CALL  PLOT  (LAMBDA, 
GO  TO  0900 
4200   CONTINUE 
STOP 
END 


)*60.+.5 

NMINX,X,NDEGY,NMINY,Y 

,«-'  ,12, Al  ,'  ,'  ,  13,'-'  ,12, Al, '     IS    THE 

LUTION. ' ) 

•"OBJECTIVE    FUNCTION    =«,F8.2,«     DEGREES', 


NUMBER 

13, «     DEGREES    OF     FREEDOM') 


OLUTICN    SIGNAL 
ZIMUTH       AZIMUTH 
X,'     NR    SITE    NAME 
P  H       N  M       A  Z  •  ) 

ER 


SITE',/, 30X, 

RANGE     PLOT', 
BETA       -       THETA     « 


3800 


3900 


6 
4003 


DDEG 

,PSI,SITE(  1,K),SITE(4,K),SITEAZ, 


) 

K)  )*DEGRAD 
LEB)*SIN(ARCC) ) 
RCB)*CGTAN(ARCC) ) 

TAT=180.-STAT 
.  )    GO    TO    3  800 
EA*RADDEG 
AZTHET=AZTHET+360. 

6.)     STAAZ(K)=STAAZ(K)-360. 

EA*RADDEG 

)    AZTHET=AZTHET-360. 

)     STAAZ(K) =STAAZ(K)+360. 

2 

E(8,K),SITE(9,K),SITE(10,K),BETA(K), 

MSITE,AZTHET,NMTHET,STAAZ(K) 

3A4,F8.3,F8.3,F7.3,I6,F6.1,I4,F6.1) 


EAN  ARC  STATISTIC  =«,F5.1,«  DEGREES. 


, NUMBER 
-STAAZ(J) ) 

ACUTE=ACUTE-180, 
ACUTE=180. -ACUTE 


T 

ACUTE  ANGLE  OF  INTERSECTION  =»,F5.1, 

PSI, DEGF, NUMBER) 
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C  GGL 


1000 
1100 
1200 


C  FIN 
C  (  0, 


C  FIN 
C  (-P 


1300 
1400 

1500 
160  0 
1700 


iaoo 

C  RET 

C  INT 

1900 


2000 


C  TES 
2100 


SURRO 
DFN  SF 
RFAL 
DIMFN 
CCMMO 
CCMMO 
FATA 
1RADDF 
SAVFL 
SAVFP 
I  FLAG 
BEGIN 
END=1 
J  =  l 
DFL(  J 
GO  TCI 
I  FLAG 
DEL{  J 
FDFL( 
LAMBD 
PSH  J 
D  GORR 
TWOPI ) 
IF  (L 
IF  (L 
D  CORR 
I  HALF, 
IF  (P 
IF  (P 
GGSLA 
SINLA 
CCS  PS 
S1NPS 
VARLA 

on  17 

CALL 
1 A  «  B  ,  C 
TEST  = 
IF  (A 
IF  (T 
C  =  C+T 
GO  TO 
C=C-T 
C^EGA 
FDEL( 
IF  (I 
J  =  2 
GG  TO 
CCNT  I 
A  IN  VA 
ERVAL 
IF  {F 
END=D 
PSI  (2 
LAMBD 
FDFL( 
DFL(2 
J-l 
GO  TO 
B  F  G  I N 
PSI  (1 
LAMBD 
FDFL( 
DEL(1 
J=2 
T  SIZE 
IF  (  ( 
GO  TO 


SEARCH* L  ,P  ,X,Y,SAVEF) 

SEARCH 
A  (  ?  )  ,  L 
DFL( 2) ,FD 
FS/S  ITE(1 
NAL/THETA 
/6. 283185 
29  57  8/ 


UTINE 
CTICN 
L  A  M  R  D 
SION 
M/SIT 
N/SIG 
TWOPI 
G/57. 
=  L 
=  P 
=  0 
=-1. 


)=0.61S*BEGIN+0. 382*END 

1200 
=  2 

}  =  0  .  38 2*B  EG  I  N+  0.  6  1  8*  E  ND 
J)=0. 

A( J) =L+X*DEL( J) 
)  =  P+Y 
FSPCN 


ELI  2)  ,PSI  (2) 

0,19)  .NUMBER 

(19) 

/, PI  HALF/ 1.57 0796/, PI/ 3. 14 1593/, 


L+X*DEL( 
r*DEL  (J) 

ICING    TEST 


LONGITUDE     IN    THE     INTERVAL 


AMBDA 
AMBCA 
ESPON 


(J)  . 
(J). 

D I  NG 
PI  HALF)  . 
SI  (J) 
S  I  (  J  ) 
M=CCS 
M=SIN 
I  =  COS 
I  =  S  I N 
M=LAM 
00  K 
AZCCM 
) 

A*B+C 
BStTE 
EST) 
W0P1 

1600 
WCPI 
=A*B+ 
J)  =FD 
FLAG. 


GT.TWCPI  )  LAMBDAT  J ) =LAMBDA ( J )-TWOP I 
LT.O.  )  LAMBDA( J ) =TWOP I +LAMBDA ( J ) 
TEST  LATITUDE  IN  THE  INTERVAL 


*  L   !  * 

(LAM 
(LAM 
(PSI 
(PSI 
BDA( 
=  1, 
P(CO 


F)     PSI ( J)=PI-PSI(J) 
LF)     PSI ( J)=-PI-PSI ( J) 
)) 
)  ) 


PIHALi 
-PIHAI 

BCA(  J 
BDA(  J 
(J  )  ) 
(J)  ) 

J) 

NUMB  El 

SLAM, SI NLAM,COSPSI,SINPSI,K, VARLAMt 


R 


-THETA(K) 
ST)-PI)  1 
1400, 1600 


C-THETA(K 
EL( J)+CME 
GT . 1  )  GO 


600,1600, 1300 
,1500 


GA**2 
TO  1800 


1100 
NUE 
LUES 
AND  T 
DFL(1 
EL(2) 
)  =  PSI 
A  (  2  )  = 
2)  =FD 
)  =  DEL 

2100 
=  DEL( 
)  =  P  S I 
At  1)  = 
1  )=FD 
)  =DEL 

OF  T 
END-B 
(  100 


CORRESPCN 
EST  NEXT 
)-FDEL(2) 

(1) 

L£MBCA(1 ) 
EL(1  ) 
(1  ) 


1) 
(2) 

LAMBDA(2) 
EL(2  ) 
(2) 

H^  INTERV 
FGIN) .LE. 
0,  1100), 


CING  TO  CURRENT  MINIMUM,  ADJUST 
DELTA  SYMMETRIC  WITH  CURRENT  DELTA. 
)  1900,1900,2000 


AL 

0.01)  GO  TO  2200 

J 
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2200   MI N=l 

IF  (J.F0.1)  MIN  =  2 

L=LAMBDA(M  IN) 
P=PSI (MINI 
DELTA=DEL( MI N) 
DEGF=FDEL(  MI N ) >  R ADDEG**2 
IF  (DEGF.LE.SAVEF)  GO  TO  2300 
L=SAVEL 
P=SAVEP 
RETURN 
2300   SAVEF=DEGF 
RFTURN 
FND 


C  A 
C  R 


100 
110 


120 

130 
140 
150 


160 
170 


SUBROUTI 

ZDTST  COMPU 
ESPFCT  TC  A 

DOUBLE  P 
1  DEN, COS B 

DATA  TWO 

REFLAM=R 

REFPSI=R 

TAPLAM=T 

TARPSI=T 

COSRP=DC 

SINRP=DS 

COSTP=DC 

SINTP=DS 

UC=DCOS( 

RNUM=SIN 

UG=SINTP 

IF  ( UG . G 

IF  (UG.L 

£RC=CARC 

DEN=DSIN 

IF  (DABS 

IF  (DABS 

CO SB =0-0 

GO  TO  11 
0  COSB=RNU 
0       CONTINUE 

CCSB=RNUM/CE 

IF     (COSB.GE. 

IF    (COSB.LE. 

BFTA=CARCPS( 

TEST=REFLAM- 

IF    (TEST)     12 

IF     (TEST+PI) 

IF     (TEST-PI) 

BETA=TWOPI-B 

AZ=BETA*RACD 

DI  ST=ARC*343 

GO    TO    1700 

AZ=0. 

DIST=0. 

CONTINUE 

RETURN 

END 


NE    AZDIS 
TES    AZIM 
RE  PER  EN 
RECISION 
.REFLAM, 
PI/6.283 
LAM 
PSI 
LAM 
PSI 

OS  (R FTPS 
IN(REFPS 
OS (TAR  PS 
INdARPS 
TARLAM)*.= 
TP*COSRP 
*SINRP+C 
E.  l.ODO) 
E. -l.ODO 
CSCUG) 
(ARC) 
(ARC). 
(RNUMJ 
DO 
00 
M/DEN 


T(RLAM,RPSI •TLAM,TPSI,AZ, DIST) 
UTH    AND    DISTANCE    OF    A    TARGET    WITH 
CE    POINT 

COSRP,SINRP,COSTP,SINTP,UC,RNUM,UG, 
REFPSI,TAPLAM,TARPSI ,  ARC 
18 5 /,RADDEG/ 57. 29578/, P 1/3. 141593/ 


I  ) 
I) 
I  ) 
I) 

ccc 

-co 

OST 

UG 

)  U 


S ( REFLAM }+DSIN(TARLAM)*DS I N( REFLAM) 

STP*SINRP*UC 

P*COSRP*UC 

=1.0D0 

G=-1.0D0 


LT 
•  G 


N 

l.OD 
-1.0 
CO  SB 
TARL 
00,  1 
140 
140 
ETA 
EG 
9.92 


.0.1D-4)  GO  TO  1600 
T.O.OD-10)  GO  TO  1000 


0)  C0SB=1.0D0 

DO)  COSB=-1.0D0 

) 

AM 

500, 1300 

0, 1500, 1500 

0,1500,1500 


DO 


COSPSI, SINPSI,K, 
RESPECT  TC  SITES 


SUBROUTINE  A ZCOM P ( CCSL AM , S I NL AM 

AZCOMP  COMPUTES  POINT  AZIMUTHS  WITH 

1  LAMBDA, A,B,C) 

REAL  N, LAMBDA 

COMMON/ SITES/SI  TEC  10,19)  ,  NUMBER 

CAT A  P  1/3. 14159 3/, T WOP  1/6. 283 18 5/, SCREEN/. 9999999/ 
A  =  l. 
C=0. 

UC=COSLAM*SITE ( 3 , K ) +S I NL AM*S ITE ( 2 , K ) 
N  =  SI NP S I *  S I T  E ( 6 , K ) -COS  P  S I*  S I T E ( 5 , K ) *UC 
UG=SINPSI*SITE ( 5,K)+C0SPSIvSITE(6,K)*UC 
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1000 
1100 
1200 

1300 


IF  (UG.GE.l.)  UG=SCREEN 

IF  (UG.LE.-l.)  UG=-SCREEN 

D=SIN(ARC0S(UG) ) 

C0SB=N/D 

IF  (C.nSR.GE.l.  ) 

IF  (C0SB.LE.-1. 

B=ARCOS(COSB) 

TEST=SITE( 1,K)- 

IF 

IF 


COSB=SCREEN 
)  COSB=-SCREEN 


LAMBDA 
(TEST)  1000,1300,1100 
(TFST+PI)  1200,1303,1300 
1200, 1300,1300 


IF  (TEST-PI) 

A=-l 

C=TWOPI 

RETURN 
END 


1000 
1100 

1 


SUBROUTINE  SHI FT(FP0ML , F ROMP , TOL , TCP , NRTYPE , DEGF , DI  ST ) 
THIS  SUBROUTINE  COMPUTES  AND  WRITES  SEARCH  PROGRESS 
DIMENSION  TYPES(2) 

DATA  RADDEG/57. 29578/,  FNCECC/1 . 006768/ 
DATA  TYPES/' GOLD' ,' GR AD'/ 

CALL  AZDIST(FROML,FRCMP,TCL,TOP,AZGQ,DIST) 
DEGLAM=TCL*RADDEG 

DFGPHI=ATAN( FNC  ECC-'TAN  (  TOP  ) )*RADDEG 
IF  (NRTYPE. EO. 2)  GO  TO  1000 
TYPF=TYPES(1) 
GO  TO  1100 
7YPE=TYPES(2) 
WRITE  (6,1)  TYPE, D I  ST , AZGO, DEGL AM , DEGPH3 ., DEGF 


FORMAT  (31X,A4,F7.1 
1F8.3,  «.  F=',F10.3) 
FROML=TOL 
FROMP=TOP 

RETURN 
END 


'  NM  ON  AZ«, F6.1,'  TO  '  , F8  .  3 , ' E ' t 


0900 


1000 
1100 


SUBROU 

PLCT  COMPU 

RPAL  L 

DOUBLE 

ICC SANG 
COMMON 
D I  MENS 
DI  f-'FNS 
DIMFNS 
FATA  S 
DATA  P 
DATA  D 

1RADDEG 
DATA  C 

114.7,1 
BCUND= 
PROB=. 
NSIZE= 
SIZE=N 
TOP=PS 
KCUNT= 
MARKV= 
DO  270 
ADD=0. 
IF  (MA 
DO  100 
OUT( J) 
CO  NT  IN 
GO  TO 
TAKF=I 
DIRECT 
VARPSI 
IF  (VA 
VARPSI 


TINE     PLOT(LAMBDA,PSI ,DEGF 
TFS    AND    WRITES    THE    CONFIDi 
AMBDA,L 
PRECISI ON    VARPSI, PI  HALF, PI, ARC, COSCCL, COL AT, 


NUMBER) 
iNCE    REGION 


/SIG 

ION 

ION 

ION 

YMBO 

IHAL 

EGRA 

/57. 

HI90 

6.  ,1 

CHI9 

900 

400 

SIZE 

1  +  1. 

0 

0 

0    I 


NAL/THETA(19) 

OUT(  52),  SYMBOLO) 

CHI900 (19 ) 

JSCALE (11 ) 

L/«*» , •     '  ,  'S'/ 

F/ 1.5  7  0796326  DO/, PI/ 3. 14 159  265 300/ 

0/1.74  53  2  9E-2/,TW0PI /6. 2  831 85/, 

29578/ 

0/2. 71, 4. 61, 6. 25, 7. 78, 9. 24, 10. 6,  12. 0,13. 4, 

7.3, 13.5,19. 8,21.1,22.3,23.5,24.8,26. ,27.2/ 

00 (NUMBER) 


/2 


05vSIZE*2.90  70443E-4 


=  1,41 


TO  1100 


RV.NE.2)  GO 
0  J  =  1,  51 
= SYMBOL ( 1) 
UF 
2500 


=  1. 

=TOP-TAKE*0. 0  5*SIZE*2.90  70443E-4 
RPSI.LT.PIHALF)  GO  TO  1200 
=PI-VARPSI 
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DIRECT=-1. 

ADD=PI 
GO  TO  1300 
1200   IF  { VARPSI .GT.-PIHALF)  GO  TO  1300 
VARPSI=-PI-VARPSI 
DIRECT=-1. 

add=pi 
130  0  parc=1 . 04«size*2.9  0  70  4  43e-4 

MARKH=0 

DO  2500  J  =  1,51 

IF  (MARKH.EC.2)  GO  TO  2100 

TAKE=J 

ARC=PARC-TAKE^0.04-*SIZE^2.9  0  704A3E-4 

COSC OL =  DCC S ( ARC) *DS IN (VARPSI ) 

IF  (COSCOL.GT. l.DO)  C0SC0L=1.D0 

IF  (COSCCL.LT.-l.DO)  COSCOL=-1.DO 

COLAT=DAPCCS (COS  COL) 

P=PIHALF-COLAT 

CCSANG=DCOTAN(VARPS I ) ADCQTAN ( COL  AT ) 

IF     (COSANG.GT. 1.D3)     C0SANG=1.D0 

IF     (COSANG.LT.-l.DO)    C0SANG=-1.D0 

ANGLE=DARCCS(CCSANG ) 

IF    (J.LE.26)    GO    TO     1400 

L=LAMBDA+DIRECT*ANGLE+ADD 

GO    TO    1500 
140  0       L  =  LAMBDA--DIRECT*ANGLE+ADD 
1500       CCNTTNUE 

IF     (L.GT.TWGPI)     L=L-TWOPI 

IF     (L.LT.O.)    L=TWOPI+L 

COSLAM=COSCL) 

SINLAM=SIN(L) 

cnspsi=cos(p ) 

SINPSI=SIN(P) 

DEGF=0. 

CO  2000  K  =  1, NUMBER 

CALL  A ZC CMP (COS LAM, SI NL AM, COS PS I,SINPSI,K,L,A,B,C) 

TEST=A#B+C-THETA(K) 

IF  (ABS(TEST)-PI )  19C 0, 1900, 1600 
1600   IF  (TEST)  1700,1900,1800 
1700   C=C+TWOPI 

GO  TO  1900 
1800   C^C-TW^PI 
1900   OMEGA=A*B+C-THETA(K) 
2  00  0   CFGF=DEGF+{CMEGA*RACDEG)**2 
2100   OUT ( J) =SYM30L( 1) 

MARKER=MARKH+1 

GOTO  (2300,2200,2500),  MARKER 
2200   IF  ( DEGF.GT. BOUND)  MARKH=2 
2300   IF  (DEGC  .GT.30UNC)  GO  TO  2400 

OUT (J) =SYM^0L(2 ) 

K0UNT=K0UNT+1 

MARKH=1 
.  MARKV=1 

IF  ( I.GT.l)  GO  TO  2350 

NSIZE=NSIZE+200 

GO  TO  0900 
2350   IF  (J.GT.l)  GO  TO  2400 

NSIZE=NSIZE+NSIZE 

WRITE  (6,7) 
7      FORMAT  (  ']  '  ) 

GO  TO  0900 
2400   IF  (I.NE.21)  GO  TO  2500 

IF  (J.E0.26)  OUT( J )=SYMB0L(3) 
2  50  0   CONTINUE 


2600 


IF  (MARKV 

.NE 

.1) 

GO  TO 

2600 

IF  (MARKh 

'•EC 

i.O  ) 

MARKV= 

2 

CCNTINUE 

IF  ( I .NE. 

1) 

GO 

TO  2650 

WRITE  (6, 

1) 

PRO 

B,NSIZE 

,NSIZE 

FORMAT  { ■ 

0«  , 

3JX 

,F6.3,« 

CONFIDENCE 

REGION 

1  •  NM  BY' , 

I  5, 

•  NM 

) 

WRITE  (6, 

3) 

PLOT  0N«,I5, 


95 


3 

2650 


4 
2700 


2800 


17C '  ),'—') 


^SIZE)+ 
■NM 

(J) 

•  I  •  .51A1,« I •  ) 


J=l,51) 


FORMAT  (37X 

CO NT  INUE 

VFRT=21-I 

NM  =  A3S(  0.05*  VERT' 

IF  (  I  .GT.21  )  NM=- 

WRITE  (6,4)  Nil,  (CUT 

FORMAT  ( »  ■ ,30X, 15, 

CONTINUE 

WRITE  (6,3) 

DO  2800  J  =  1,11 

ADD  =  J 

JSCALE( J) =ABS( (0.2-ADD-1.2KSIZE)  +  .5 

IF  (J.LT.6)  JSCALEl J)=-JSCALE( J) 

WRITE  (6,5)  (JSCALE(J)  ,J=1,11) 

FORMAT  (34X.11I5) 

IJAREA=0.002ASIZF**2+.5 

N  M AR  E A  =  KCUNT* IJ  A  R  E  A 

WRITE  (6,6)  NMAREA 

FORMAT  (30X,«      CONFIDENCE  REGION  AREA 

«  NAUTICAL  MILES  SQUARED.') 

RETURN 

END 


17, 
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GOLD 

4  398 

.4 

NM 

ON 

AZ 

270.0 

TO 

2  86. 

74  OE 

0.0 

F  = 

20254. 

277 

GOLD 

4675. 

,4 

NM 

ON 

AZ 

180.0 

TO 

286. 

,740E 

-77.953. 

F= 

11288. 

.250 

GOLD 

1021. 

.9 

NM 

ON 

AZ 

134.1 

TO 

16. 

323E 

-77.953. 

F= 

8562 

.371 

GOLD 

1898 

.3 

NM 

ON 

AZ 

0.0 

TO 

16. 

323E 

-46.449. 

F  = 

3712. 

.686 

GOLD 

1108. 

,5 

NM 

ON 

AZ 

99.8 

TO 

43. 

155E 

-46.449. 

F  = 

868. 

.7  00 

GOLD 

412. 

.8 

NM 

ON 

AZ 

0.0 

TO 

43. 

155E 

-39.571. 

F  = 

Ill 

.312 

GRAD 

271 

.3 

NM 

ON 

AZ 

94.0 

TO 

49. 

004E 

-39.741. 

F  = 

15. 

.526 

GRAD 

84. 

.6 

NM 

ON 

AZ 

53.8 

TO 

50. 

461E 

-38.898 . 

F  = 

0 

.158 

GRAD 

10. 

.2 

NM 

ON 

AZ 

76.8 

TO 

50. 

673E 

-33.859. 

F  = 

0. 

.050 

GRAD 

0. 

.1 

NM 

ON 

AZ 

52.7 

TO 

50. 

675E 

-38.858. 

F= 

0. 

.050 

50-40F,     38-51S    IS    THE    LEAST    SQUARES     SOLUTION. 
OBJECTIVE     FUNCTION    =  0.05    DEGREES    SQUARED. 

3    DEGREES    OF    FREEDOM 


SOLUTION 

SIGNAL 

SITE 

ID 

INTERCEPT 

AZIMUTH 

AZIMUTH 

RANGE 

NR 

SITE  NAME 

BETA   - 

THETA  = 

OMEGA 

NM 

9 

HONOLULU 

228.098 

228.241 

-0. 142 

9001 

19 

WICHITA 

102.627 

102.797 

-0.170 

9301 

10 

LOS  ANGELES 

121.457 

121.424 

0.033 

10199 

PLOT  THETA  VALS 
NORM   NM   AZ 

27.1  4    117.1 

187.3         4    277.3 
335.1         0    245.1 

MFAN  ARC  STATISTIC  =  21,8  DEGREES. 

MEAN  ACUTE  ANGLE  GF  INTERSECTION  =  34.7  DEGREES. 

0.900  CONFIDENCE  REGION  PLOT  ON   400  NM  BY   400  NM  AREA 


200 

190 

180 

170 

160 

150 

140 

130 

120 

110 

100 

90 

80 

70 

60 

50 

40 

30 

20 

10 

0 

-10 

-20 

-30 

-40 

-50 

-60 

-70 

-80 

-90 

■100 

-110 

■120 

■130 

-140 

■150 

-160 

-170 

■180 

-190 

■200 


& -i sfe, & sfe <*: .'<'J  :>- t^$$$ ?  ft$'$  t  $  ft#&*  >$  $ 
*s  r*';*:;"  r  ■?'  -t'.'TV.*  :'•:r.-.•^," :;":-'  :'.  :.'  :  -  f  ;'  :;  :"•■  •$  :'  ;' 
^*  •(:"■•.;     '   -■:[    $%%&$:&&&.#&.&'#£:&&$ 
:V-:^-''  :    :.      ■-  >i  :'•  s\  £.  »;  >  .'.  :'  .-  i  i.  t:  £  &'Jz&i*. 
>"*>%' :*:-::^i'::*;YX'.:- ■  :'■  ■:;:    :'  ^:k?>  sesjssj  t  ■  *.*•;- 

^•^  ::-^  'ifij^:-  ^;:/::  '•-.\  i-.:'r  &&  sfc&jjsjj  :* 
*£* 

:H ■*&*:%;& *&&**  ft-  *  fc;X* 


<  ^ 


.  >  .  j. ..  .  - 


*  r!   | 

•>'  .-    :: 

r*  r*  ■ : 

.•-  -w;.- 

-•'  *  $f 

*#* 


*  i;  ft  •' 

v  :;-  sr  :r. 
*y.  If  * 


r  *  fe 

-.-  .'■    • 

s  *;*£ 

S;  A  « 

:rf-  y 

>■?  rfr  fe 
sfe&jfe 

s*  .<  6 


t    y  *  A 

i*  j»  :;-;  ?, 

■fS  A  ;<^  ft' 
$i  F  .t  A' 
*6  *  ":  V- 

*  •?  i  ^» 


i;  i<.  v. '.':  f-  J 
:■■  ■•*::*::■,  *  J 
j  >4:^  I 
*#- 1**  I 
rV^ryr*:;:! 

:'-«■  *$*! 
:;.:';.  si  J 
^ fet** I 


1 .  rf- : *•  •---  i[ .  % :"  *  £  &  ft  ?•..  *  i;  t  ft  sf.  t ■.  is  A  *  ^  -^  ^s  •> :*-  S . s% : t  ,'s  *  *  >:  A  :»  *  $  *.  *  -   :'<■  4'  A  *  sfe : 


5  5^,  A: 


--.:;  ■,•:-: 


:■  r/A^i^;- 


;<--V^  ;..- ■'  ••'.'■.'-: ;<t:*->i  ^  i. -5r  i' -J- ?;•  <-^^Sf 


•t:'. 


^•':(t'  :fejfv!%; 


*sfe^r%±*A» 


>  r;  r>  A  i  A  $  *  :*'. 


S  .   ;•  *  ;   i.-  $  r"  v"  fe  *  *  ^'^  ^  ^  •*  >  *  ^«  s 


A*'  ; 


:*: 


'^y-*:f;A: 


^?;> 


.r<f'->  * 


$^!^.^:AAsSe*  I  •*. 


±& 


\t  i  S._  ^  &  $  A  !*  5fe  *  'if.  $  ^t  Jj£  SfK 

$"< s«-  *•  *;> ^; re ^ :t- ^  ^? "t ~$. '  •  * *■* !* ' 
*•  r»  :*  *;  •*? ;  i  "ii>  K-  -;*.  tit  A  f' .  $■  4- 
^;  f"  v.^!  ■&  '4: . ?;i;  ^:*  «r  fe  if 


ft  *c  ft  *•  fe *  H  &  -- : 


\  T 


^  ?-■  ;r'»  j ; 


•5*.^i-*fi;*  t?5" 


ft:  ^.?  S-  ^'  'h  -^  :••  &  a  ft  $  :*=  .v  s5  ri 


T'<^.T 


: ^rf:fe;i>  rfe-ft:'- rfr&: 


>;.  l£  $  4;  st  f>>.  r*r  *  <;  ^  -^t  j£  "i  •*;  i;;^r  :*  ar.  r^?4 : 


•4  jfc.  rt  «s  ft  4.  at :?.  ^  *  t  ii  ar  J 

^  fc  $  rf-  ^    15?  ^ r-  ^"  «  *    J 


,<:u:^^.--:S:r 


:  A  ft  ^i :  s  w  <-  ^-  J'-  *5f  :*  ;•%  #:  is  Tf; 


A  ft  °«s*S'  *^  *  4s*$A  r'-  ^-  K?  AAAv^  •*  $  iteA^^^rjfc*  rf;^:*^  ^  A  £  *  ±  ,->  *:s 


-&: 


rsfesSssJi: 


■•..  ri  $  }i:  5  r-5  *  A  j.:  ¥  iy  r%  :&  >.  >;  ^  *;  A:  ^  rfe  *  :i  if  £■  &  *  st"  H': 


>«<^X  ^    ^  f>  A  ^  ^    J 


J 


^■:$ft  »-:'ct--:!TSpr:;^r-iV^-s;';j'::J5:'-.^^^-^  ^  ft  si;  :*  >^  <s  ti  &  rl-  &  A ^^  ife  A  *.  ^  *  *  *  tS-  *  (t*$*  ?  ■£  ££:;**■  1 
;f  fc  r*:  rk  ^  r*^  4  i  r$  r>  *  jft  r-feijert  S"  ^  *  sf :?  :•  r*-  tf:  A  ri  »  S;  A  $  *  :'-  & : ';  v  :y  :•?  ^t  jL  S  i  »  v  $  tt'^>  *  t-  *  * :>  I 


-200  -160  -120   -80   -40     0    40    80 

CONFIDENCE  REGION  AREA  =    6080  NAUTICAL 


120   160   200 
MILES  SQUARED. 


97 


GOLD 

6104.5 

NM 

ON 

AZ 

270.0 

TO 

258. 

323E 

0.0 

F= 

86.016 

GOLD 

276.9 

NM 

ON 

AZ 

180.0 

TO 

258. 

323E 

-4.643. 

F= 

79.361 

GRAD 

2293.2 

NM 

ON 

AZ 

167.1 

TO 

268. 

963E 

-41 .796. 

F= 

10.165 

GPAD 

1199.7 

NM 

ON 

AZ 

160.2 

TO 

282. 

300E 

-59.978. 

F  = 

0.043 

GRAD 

77.3 

NM 

ON 

AZ 

151.7 

TO 

283. 

557E 

-61.102. 

F= 

0.013 

GRAD 

0.9 

NM 

ON 

AZ 

97.3 

TO 

283. 

589E 

-61.104  . 

F= 

0.013 

GRAD 

0.0 

NM 

ON 

AZ 

0.0 

TO 

283. 

590E 

-61.104. 

F  = 

0.013 

76-25W,  61-  6S  IS  THE  LEAST  SQUARES  SOLUTION. 
OBJECTIVE  FUNCTION  =  0.01  DEGREES  SQUARED. 
3  DEGREES  OF  FREEDOM 


ID 
NR 

7 
17 

8 


INTERCEPT 
SITE  NAME 
EL  PASO 
RENO 
HELENA 


SOLUTION 

AZIMUTH 

BETA   - 

165.841 

159.655 

162.349 


SIGNAL 
AZI MUTH 
THETA  = 
165.922 
159.655 
162.270 


OMEGA 
-0.081 
0.000 
0.080 


SITE 

RANGE 

NM 

5747 

6383 

6673 


PLOT  THETA  VALS 
NORM   NM   AZ 
244.6    4  334.6 

56.4  0  326.4 

64.5  4  334.5 


MEAN  ARC  STATISTIC  =  75.6  DEGREES. 
MEAN  ACUTE  ANGLE  OF  INTERSECTION  = 


5.5  DEGREES. 


0.900  CONFIDENCE  REGION  PLOT  ON  2400  NM  BY  2400  NM  AREA 
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480    I 

420    1 

360    I 

300    I 

240     I 

180    I 

120    I 

60    I 

0    I 

-60    I 

-120    I 

-180    I 

-240     I 
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-1200  -960  -720  -480  -240     0   240   480 
CONFIDENCE  REGION  AREA  =  311040  NAUTICAL 


720   960  1200 

MILES  SQUARED. 


98 


GOLD 

7810.8 

NM 

ON 

AZ 

270.0 

TO 

229.901E 

0.0   . 

F= 

663.993 

GOLD 

995.3 

NM 

ON 

AZ 

180.0 

TO 

229.901E 

-16.684. 

F= 

13.884 

GOLD 

97  .0 

NM 

ON 

AZ 

90.2 

TO 

231. 588F 

-16.684. 

F  = 

3.408 

GOLD 

3.8 

NM 

ON 

AZ 

0.0 

TO 

231.588E 

-16.620. 

F= 

3.381 

GRAD 

10.1 

NM 

ON 

AZ 

67.7 

TO 

231.750E 

-16.556. 

F  = 

3.289 

GRAD 

0  .0 

NM 

ON 

AZ 

0.0 

TO 

231.750E 

-16.556. 

F  = 

3.289 

128-15W,  16-33S  IS  THE  LEAST  SQUARES  SOLUTION. 
OBJECTIVE  FUNCTION  =  3.29  DEGREES  SQUARED. 
4  DEGREES  OF  FRFEDOM 


SOLUTION 

SIGNAL 

SITE 

ID 

INTERCEPT 

AZIMUTH 

AZIMUTH 

RANGE 

PLOT  ' 

FHETA  VALS 

NR 

SITE  NAME 

BETA   - 

THETA  = 

OMEGA 

NM 

NORM 

NM   AZ 

14 

NASHVILLE 

224.264 

224.377 

-0. 113 

3929 

306.1 

6   36.1 

2 

ATTU 

124.892 

126.340 

-1.44  8 

5135 

238.9 

86  328.9 

9 

HONOLULU 

140.026 

138.983 

1  .043 

2854 

50.6 

46  320.6 

16 

NORFOLK 

232.077 

231.774 

0.303 

4406 

130.0 

17   40.0 

MEAN  ARC  STATISTIC  =  68.0  DEGREES. 

MFAN  ACUTE  ANGLE  OF  INTERSECTION  =  50.9  DEGREES. 

0.900  CONFIDENCE  REGION  PLOT  ON   400  NM  BY   400  NM  AREA 
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20880  NAUTICAL  MILES  SQUARED. 


99 


GCLD    4399.1    NM    ON    AZ       90.0    TO       73.272E 


0.0 


GOLD 

1765.6 

NM 

ON 

AZ 

0.0 

TO 

73.272E 

29.574. 

F 

GOLD 

123.9 

NM 

ON 

AZ 

89.4 

TO 

75.642E 

29.574. 

F 

GOLD 

71  .1 

NM 

ON 

AZ 

0.0 

TO 

75.642E 

30.763. 

F 

GRAD 

14.2 

NM 

ON 

AZ 

55.7 

TO 

75.869E 

30.897. 

F 

GRAD 

0.2 

NM 

ON 

AZ 

40.6 

TO 

75.871E 

30.899. 

F 

GRAD 

0.0 

NM 

ON 

AZ 

0.0 

TO 

75.871E 

30.899. 

F 

75-52F,  30-54N  IS  THE 
OBJECTIVE  FUNCTION  = 
5  DEGREES  OF  FREEDOM 


LEAST  SQUARES  SOLUTION. 
3.30  DEGREES  SQUARED. 


767.3M 
26.527 
4.538 
3.487 
3.303 
3.303 
3.303 


SOLUTION 

SIGNAL 

SITE 

ID 

INTERCEPT 

AZIMUTH 

AZIMUTH 

RANGE 

PLOT  ' 

rHETA 

VALS 

NR 

SITE  NAME 

BETA   - 

THETA  = 

OMEGA 

NM 

NORM 

NM 

AZ 

17 

RENO 

345.832 

345.356 

0.476 

6511 

102.9 

27 

12.9 

5 

DENVER 

359.185 

358.897 

0.288 

6589 

90.8 

16 

0.8 

12 

MIDWAY 

304.587 

304.956 

-0.369 

5329 

327.7 

22 

57.7 

19 

WICHITA 

6.275 

7.726 

-1.451 

6693 

263.  7 

81 

353.7 

7 

EL  PASO 

357.715 

356.849 

0  .866 

7064 

92.7 

46 

2.7 

MEAN  ARC  STATISTIC  =  72.3  DEGREES. 

MEAN  ACUTE  ANGLE  OF  INTERSECTION  =  28.0  DEGREES. 

0.900  CONFIDENCE  REGION  PLOT  ON   400  NM  BY   400  NM  AREA 
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.ONFIDENCE  REGION  AREA  =   36480  NAUT 


80   120   160   200 
ICAL  MILES  SQUARED. 


100 


GGLD  2476.5  NM  ON  AZ   90.0  TO   41.249E 


0.0 


.  F  = 


GOLD 

561.3 

NM 

ON 

AZ 

0. 

0 

TO 

41. 249E 

9.411. 

F 

GOLD 

2  07.1 

NM 

ON 

AZ 

270. 

3 

TO 

37.754E 

9.411. 

F 

GOLD 

213.3 

NM 

ON 

AZ 

0. 

0 

TO 

37.754E 

12.985. 

F 

GRAD 

277.8 

NM 

ON 

AZ 

318. 

0 

TO 

34.526E 

16.421. 

F 

GRAD 

4.4 

NM 

ON 

AZ 

325. 

8 

TO 

34. 48 3 E 

16.482. 

F 

GRAD 

0.  1 

NM 

ON 

AZ 

325. 

7 

TO 

34.482E 

16.484. 

F 

GRAD 


0.0  NM  ON  AZ 


0.0  TO   34.482E   16.484.  F= 


258 

.677 

94. 

.282 

42 

.317 

20 

.378 

0, 

.264 

0 

.259 

0 

.259 

0 

.259 

34-29E,  16-29N  IS  THE  LEAST  SQUARES  SOLUTION. 
OBJECTIVE  FUNCTION  =  0.26  DEGREES  SQUARED. 
5  DEGREES  OF  FREEDOM 


SOLUTION 

SIGNAL 

SITE 

ID 

INTERCEPT 

AZIMUTH 

AZIMUTH 

RANGE 

PLOT  ■ 

rHETA 

NR 

SITE  NAME 

BETA   - 

THETA  = 

OMEGA 

NM 

NORM 

NM 

AZ 

8 

HELENA 

34.383 

33.967 

0.417 

6624 

66.2 

23 

336.2 

11 

MIAMI 

63.878 

64.000 

-0.122 

6235 

212.4 

7 

302.4 

14 

NASHVI  LLE 

57.577 

57.539 

0.038 

6229 

44.6 

2 

314.6 

17 

RENO 

28.527 

28.765 

-0.238 

7165 

247.2 

12 

3  37.2 

1 

ANCHORAGE 

355.710 

355.825 

-0.114 

6153 

272.1 

6 

2.1 

MEAN  ARC  STATISTIC  =  72.1  DEGREES. 

MEAN  ACUTE  ANGLE  OF  INTERSECTION  =  28.4  DEGREES. 

0.900  CONFIDENCE  REGION  PLOT  ON   600  NM  BY   600  NM  AREA 


300  I 

285  I 

270  I 

255  I 

240  I 

225  I 

210  I 

195  I 

180  I 

165  I 

150  I 

135  I 

120  I 

105  I 

90  I 

75  I 

60  I 

45  I 

30  I 

15  1 

0  I 

-15  I 

-30  I 

-45  I 

-60  I 

-75  I 

-90  I 

105  I 

120  I 

135  I 

150  I 

165  I 

180  I 

195  I 

210  I 

225  I 

240  I 

255  I 

270  I 

285  I 

300  I 
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57600  NAUTICAL  MILES  SQUARED. 


101 


GOLD 

1690  .8 

NM 

ON 

AZ 

270.0 

TO 

331.837E 

0.0 

F  = 

102.897 

GOLD 

141.  5 

NM 

ON 

AZ 

180.0 

TO 

331.837E 

-2.373. 

F  = 

83.012 

GRAD 

629.1 

NM 

ON 

AZ 

135.7 

TO 

339.241E 

-9. 884. 

F  = 

3.8  63 

GRAD 

2.2 

NM 

ON 

AZ 

51.3 

TO 

339.270E 

-9.861. 

F  = 

3.846 

GRAD 

0.  1 

NM 

ON 

AZ 

48.4 

TO 

339.271E 

-9.860. 

F  = 

3.845 

20-44W.   9-52S  IS  THE 
OBJECTIVE  FUNCTION  = 
6  DEGREES  OF  FREEDOM 


LEAST  SQUARES  SOLUTION. 
3.85  DEGREES  SQUARED. 


ID  INTERCEPT 
NR  SITE  NAME 


SOLUTION  SIGNAL 
AZIMUTH   AZIMUTH 


SITE 

RANGE    PLOT    THETA    VALS 


BETA 


THETA    =    OMEGA       NM 


NORM   NM   AZ 


3 

BOSTON 

125.857 

125  .554 

0.303 

4166 

52.3 

17 

322.3 

12 

MIDWAY 

54.831 

54.455 

0.376 

9092 

43.3 

10 

313.3 

2 

ATTU 

19.885 

19.173 

0.  712 

8140 

78.4 

29 

348.4 

11 

MIAMI 

113.543 

114.172 

-0 .629 

4071 

213.2 

34 

303.2 

1 

ANCHORAGE 

58.837 

60.428 

-1.591 

7009 

244.4 

85 

334.4 

10 

LOS  ANGELES 

94.063 

93.642 

0.421 

6103 

32.9 

24 

302.9 

MFAN  ARC  STATISTIC  =  58.6  DEGREES. 

MEAN  ACUTE  ANGLE  OF  INTERSECTION  =  22.0  DEGREES. 

0.900  CONFIDENCE  REGION  PLOT  ON   400  NM  BY   400  NM  AREA 
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-200  -160  -120   -80   -40 
CONFIDENCE  REGION  AREA  = 


0    40    80   120   160   200 
26720  NAUTICAL  MILES  SQUARED. 


102 


GGLD 

7158.7 

NM 

ON 

AZ 

90.0 

TO 

119.236E 

0.0   . 

F  = 

318.482 

GOLD 

251.0 

NM 

ON 

AZ 

0.0 

TO 

119.236E 

4.209. 

F  = 

261.669 

GRAD 

1861  .3 

NM 

ON 

AZ 

29.7 

TO 

136.486E 

30.744. 

F  = 

34. 133 

GPAD 

525.6 

NM 

ON 

AZ 

47.2 

TO 

144.45  IE 

36.464. 

F  = 

8.503 

GRAD 

104.9 

NM 

ON 

AZ 

34.4 

TO 

145.701E 

37.901. 

F  = 

8.042 

GRAD 

5.3 

NM 

ON 

AZ 

222.2 

TO 

145.627E 

37.836. 

F  = 

8.041 

GPAD 

0.6 

NM 

ON 

AZ 

222.0 

TO 

145.618E 

37.829. 

F  = 

8.041 

GRAD 

0.4 

NM 

ON 

AZ 

221.9 

TO 

145. 612E 

37.824. 

F  = 

8.041 

GRAD 

0.  1 

NM 

ON 

AZ 

221.8 

TO 

145.611E 

37.822. 

F  = 

8.041 

GRAD 

0.0 

NM 

ON 

AZ 

0.0 

TO 

145.611E 

37.822. 

F  = 

8.041 

145-37E.  37-49N  IS  THE 
OBJECTIVE  FUNCTION  = 
6  DEGREES  OF  FREEDOM 


LEAST  SQUARES  SOLUTION. 
8.04  DEGREES  SQUARED. 


ID  INTERCEPT 
NR  SITE  NAME 
16  NORFOLK 

4  CORPUS  CHRI 
11  MIAMI 
18  SEATTLE 

6  DETROIT 


SOLUTION 
AZI MUTH 
BETA   - 

327.990 
;  315.081 
324.266 
298.750 
323.507 


SIGNAL 

AZIMUTH 

THETA  = 

328.189 

314.619 

322.103 

300.372 

324.197 


OMEGA 
-0.199 
0.46  2 
2.  164 
-1.622 
-0.690 


SITE 
RANGE 
NM 
5685 
5524 
6215 
3873 


PLOT  THETA  VALS 

NORM  NM  AZ 

301.6  11  31, 

142.2  27  52 

132.2  126  42, 

319.2  87  49, 

303.9  41  33 


13  MINNEAPOLIS   316.823  316.840  -0.016   4923  307.8 

MEAN  ARC  STATISTIC  =  80.8  DEGREES. 

MFAN  ACUTE  ANGLE  OF  INTERSECTION  =  10.2  DEGREES. 


0   37.8 


0.900  CONFIDENCE  REGION  PLOT  ON   600  NM  BY   600  NM  AREA 


300 

285 

270 

255 

240 

225 

210 

195 

180 

165 

150 

135 

120 

105 

90 

75 

60 

45 

30 

15 

0 

-15 

-30 

-45 

-60 

-75 

-90 

-105 

120 

-135 

-150 

•165 

-180 

-195 

-210 

-225 

-240 

-255 

-270 

-285 

-300 
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-300  -240  -180  -120   -60     0    60   120 
CONFIDENCE  REGION  AREA  =   37080  NAUTICAL 


180   240   300 
MILES  SQUARED. 
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GOLD 

2341 

.5 

NM 

ON 

AZ 

90.0 

TO 

39. 

OOOE 

0.0 

F^ 

5868.273 

GOLD 

575. 

.9 

NM 

ON 

AZ 

18  0.0 

TO 

39. 

OOOE 

-9.656  . 

F  = 

5288.270 

GOLD 

486 

.3 

NM 

ON 

AZ 

90.7 

TO 

47. 

215E 

-9.656. 

F  = 

4666.781 

GOLD 

531 

.3 

NM 

ON 

AZ 

180.0 

TO 

47, 

215E 

-18.558. 

F  = 

3902.349 

GOLD 

406 

.6 

NM 

ON 

AZ 

91.  1 

TO 

54. 

355E 

-18.558. 

F  = 

3082.319 

GOLD 

425 

.4 

NM 

ON 

AZ 

180.0 

TO 

54. 

355E 

-25.695. 

F  = 

2172.056 

GRAD 

577 

.5 

NM 

ON 

AZ 

138.9 

TO 

61. 

8  49E 

-32.774. 

F  = 

215.962 

GRAD 

95. 

.2 

NM 

ON 

AZ 

36.  1 

TO 

62. 

94  2  E 

-31.483. 

F  = 

4.374 

GRAD 

10 

.9 

NM 

ON 

AZ 

99.8 

TO 

63. 

151E 

-31.514. 

F  = 

2.152 

GRAD 

0. 

.  1 

NM 

ON 

AZ 

25.4 

TO 

63. 

152E 

-31.512. 

F  = 

2.152 

GRAD 

0 

.0 

NM 

ON 

AZ 

0.0 

TO 

63. 

152E 

-31.512. 

F  = 

2.152 

63-  9E,  31-31S  IS  THE 

OBJECTIVE  FUNCTION  = 
7  DEGREES  OF  FREEDOM 


LEAST  SQUARES  SOLUTION. 
2. 15  DEGREES  SQUARED. 


ID  INTERCEPT 
NR  SITE  NAME 


SOLUTION  SIGNAL 
AZIMUTH   AZIMUTH 


SITE 

RANGE  PLOT  THETA  VALS 


BETA 


THETA  =  OMEGA   NM 


NORM   NM   AZ 


5 

DENVER 

53.228 

53 

,230 

-0.002 

10049 

223.7 

0 

313.7 

15 

NEW  ORLEANS 

100.649 

101. 

.286 

-0.637 

942  5 

182.5 

14 

272.5 

6 

DETROIT 

79.209 

78 

.237 

0.972 

9070 

32.3 

28 

302.3 

16 

NORFOLK 

90.386 

89. 

,567 

0.820 

8779 

24.3 

27 

294.3 

7 

EL  PASO 

91.177 

91 

.369 

-0.  192 

10276 

184.  1 

1 

274.1 

17 

RENO 

342.343 

342 

,395 

-0.051 

10306 

285.9 

0 

15.9 

8 

HELENA 

15.592 

15. 

,892 

-0.300 

9875 

257.2 

4 

347.2 

MEAN  ARC  STATISTIC  =  18.7  DEGREES. 

MEAN  ACUTE  ANGLE  OF  INTERSECTION  =  42.9  DEGREES. 

0.900  CONFIDENCE  REGION  PLOT  ON   400  NM  BY   400  NM  AREA 
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0    40    80   120   160   200 
1680  NAUTICAL  MILES  SQUARED. 
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GOLD 

5019.9 

NM 

ON 

AZ 

270.0 

TO 

276.387E 

0.0   . 

F  = 

348.028 

GOLD 

819  .4 

NM 

ON 

AZ 

100.0 

TO 

276.387E 

-13.737. 

F  = 

15.774 

GOLD 

73.1 

NM 

ON 

AZ 

90.1 

TO 

277.641E 

-13.737. 

F= 

5.151 

GOLD 

71.  1 

NM 

ON 

AZ 

180.0 

TO 

277.641F 

-14.929. 

F  = 

2.936 

GRAD 

1  .4 

NM 

ON 

AZ 

267.9 

TO 

277.616E 

-14.93  0. 

F  = 

2.927 

GRAD 

0.1 

NM 

ON 

AZ 

270.0 

TO 

277.615E 

-14.930. 

F  = 

2.926 

82-23W,     14-56S     IS    THE    LEAST    SQUARES    SOLUTION. 
OBJECTIVE    FUNCTION    =  2.93    DEGREES    SQUARED. 

8    DEGREES    OF     FREEDOM 


SOLUTION 

SIGNAL 

SITE 

ID 

INTERCEPT 

AZIMUTH 

A  Z I MUT  H 

RANGE 

PLOT  Th 

NR 

SITE  NAME 

BETA   - 

THETA  = 

OMEGA 

NM 

NORM 

NM 

AZ 

9 

HONOLULU 

109.222 

109.844 

-0.622 

4943 

204.5 

36 

294.5 

17 

RENO 

139.373 

138.901 

0.472 

3868 

58.4 

25 

328.4 

6 

DETROIT 

179.229 

179. 163 

0.066 

3421 

89.4 

3 

359.4 

14 

NASHVI LLE 

174.501 

173.573 

0.928 

3060 

84.8 

43 

354.8 

3 

BOSTON 

192.934 

193.473 

-0.539 

3479 

280.2 

27 

10.2 

11 

MIAMI 

183.252 

184.030 

-0.778 

2432 

273.6 

30 

3.6 

19 

WICHITA 

162.091 

161  .478 

0.613 

3253 

75.0 

29 

345.0 

8 

HELFNA 

148.667 

149.088 

-0.421 

4010 

248.4 

23 

338.4 

MEAN    ARC    STATISTIC    =    59.3    DEGREES. 

MEAN    ACUTE    ANGLE    OF    INTERSECTION    =    27.8    DEGREES. 

0.900    CONFIDENCE    REGION    PLOT    ON       400    NM    BY       400    NM    AREA 
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CONFIDENCE  REGION  AREA  = 


0    40    80 
33040  NAUTICAL 


120   160   200 

MILES  SQUARED. 
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GOLD  9485.0  NM  ON  AZ 

GOLD  74.4  NM  ON  AZ 

GP.AD  1077.6  NM  ON  AZ 

GRAD  13.2  NM  ON  AZ 

GPAD  0.0  NM  ON  AZ 

172-18F,  12-36N  I S  THE 
OBJECTIVE  FUNCTION  = 
9  DEGREES  OF  FREEDOM 


90.0  TO  157.984E 

0.0  TO  157.984F 

51.3  TO  172.232E 

15.4  TO  172.292F 
0.0  TO  172.292E 


LEAST  SQUARES  SOLUTION. 
5.86  DEGREES  SQUARED. 


0.0 

F  = 

215.253 

1.248. 

F  = 

207.475 

12.394. 

F= 

6.162 

12.607. 

F  = 

5.864 

12.607. 

F  = 

5.8  64 

SOLUTION 

SIGNAL 

SITE 

ID 

INTERCEPT 

AZIMUTH 

AZIMUTH 

RANGE 

PLOT  TF 

VALS 

NR 

SITF  NAME 

BETA   - 

THETA  = 

OMEGA 

NM 

NORM 

NM 

AZ 

9 

HONOLULU 

257.935 

258.530 

-0. 594 

1788 

339.6 

17 

69.6 

10 

LOS  ANJGELES 

269.320 

268.506 

0.814 

3968 

147.9 

44 

57.9 

11 

MIAMI 

289.107 

288.596 

0.511 

5995 

150.9 

30 

60.9 

12 

MIDWAY 

213.947 

213. 857 

0.090 

1098 

120.2 

1 

30.2 

13 

MINNEAPOLIS 

282.050 

282.304 

-0.2  54 

5064 

315.3 

15 

45.3 

14 

NASHVI LLE 

286.504 

286.771 

-0.267 

5480 

322.6 

16 

52.6 

15 

NEW  ORLEANS 

284.612 

282.885 

1.727 

5418 

149.3 

103 

59.3 

16 

NORFOLK 

293.348 

293.874 

-0.526 

5903 

317.3 

31 

47.3 

17 

RENO 

265.894 

266.982 

-1. 088 

3907 

322.7 

59 

52.7 

MEAN  ARC  STATISTIC  =  67.1  DEGREES. 

MEAN  ACUTE  ANGLE  OF  INTERSECTION  =  13.0  DEGREES. 

0.900  CONFIDENCE  REGION  PLOT  ON   400  NM  BY   400  NM  AREA 
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18080  NAUTICAL  MILES  SQUARED. 
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GO  I  D 

202. 

,  5 

NM 

ON 

AZ 

90.0 

TO 

3. 

373E 

0.0 

F  = 

10268.355 

GOLD 

1955 

.9 

NM 

ON 

AZ 

0.0 

TO 

3. 

373E 

32.769. 

F  = 

5296.391 

GOLD 

1314, 

.5 

NM 

ON 

AZ 

277.1 

TO 

337. 

319E 

32.769. 

F= 

3449.257 

GOLD 

585. 

.4 

NM 

ON 

AZ 

0.0 

TO 

337. 

319E 

42.536. 

p  = 

2502  .434 

GOLD 

2030 

.4 

NM 

ON 

AZ 

286.1 

TO 

290. 

9  72E 

42.536. 

F  = 

119.602 

GOLD 

56. 

.  1 

NM 

ON 

AZ 

0.0 

TO 

290. 

972  E 

43.471  * 

F  = 

51.410 

GRAD 

107. 

.2 

NM 

ON 

AZ 

69.9 

TO 

293. 

299E 

44.062. 

F  = 

10.708 

GRAD 

4. 

.2 

NM 

ON 

AZ 

216.4 

TO 

293. 

241E 

44.006. 

F  = 

10.553 

GRAD 

0. 

.  1 

NM 

ON 

AZ 

215.9 

TO 

293. 

240E 

44.005. 

F  = 

10.553 

66-46W,  44-  ON  IS  THE 
OBJECTIVE  FUNCTION  = 
9  DEGREES  OF  FREEDOM 


LEAST  SQUARES  SOLUTION. 
10.55  DEGREES  SQUARED. 


SOLUTION 

SIGNAL 

SITE 

ID 

INTERCEPT 

AZIMUTH 

AZIMUTH 

RANGE 

PLOT    Tt- 

NR 

SITF    NAME 

BETA       - 

THETA    = 

OMEGA 

NM 

NORM 

NM 

AZ 

6 

DETROIT 

76.495 

77.395 

-0.900 

721 

178.5 

11 

268.5 

12 

MIDWAY 

42.761 

41.147 

1.614 

5054 

33.7 

96 

303.7 

18 

SEATTLE 

74.220 

72.903 

1.317 

2294 

24.  5 

48 

294.5 

5 

DENVER 

68.794 

67.963 

0.832 

1719 

4.  1 

23 

274.1 

11 

MIAMI 

27.603 

26.907 

0.696 

1274 

304.7 

15 

214.7 

17 

RENO 

65.683 

66.781 

-1.098 

2357 

193.2 

41 

283.2 

4 

CORPUS    CHRIS       48.511 

47.035 

1.476 

1765 

335.6 

43 

245.6 

10 

LOS    ANGELES 

60.122 

60.990 

-0. 868 

2440 

184.6 

33 

274.6 

16 

NORFOLK 

50.167 

50.470 

-0.303 

534 

146.7 

2 

236.7 

MEAN  ARC  STATISTIC  =  33.6  DEGREES. 

MEAN  ACUTE  ANGLE  OF  INTERSECTION  =  33.9  DEGREES. 

0.900  CONFIDENCE  REGION  PLOT  ON   400  NM  BY   400  NM  AREA 
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